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Chapter I 
INTRODUCTION 

1.1 SURVEY OF RECENT RESEARCH 

In recent years the critical nature of the photochemical 
interactions involving chlorine, nitrogen, oxygen, and hydro- 
gen compounds occurring in the atmosphere has become apparent. 
The danger of anthropogenic reduction of the protective 
ozone layer of the stratosphere (1) and the possible global 
heating due to increases in CC >2 from deforestation and fossil 
fuels and methane from agricultural sources (2) are also of 
concern. The uncertainty surrounding the importance of each 
of these effects has stirred much interest in understanding 
the structure of the atmosphere. The study of the atmos- 
phere has recently been recognized as one of the most 
important fields of scientific inquiry both in its own right 
and because of its importance to the survival of mankind (3) . 

In order to understand atmospheric processes it is 
necessary to know the distribution of the various atmospheric 
molecules with altitude. The nature and influence of the 
pollutants in the atmosphere can only be estimated if these 
distributions are known. Both minor (parts per million) and 
trace (parts per billion) constituent distributions have 
been determined by in situ experiments which sample the local 
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environment and remote-sensing techniques which allow 
simultaneous concentration measurements over wide altitude, 
latitude, and longitude ranges of the atmosphere. 

A brief summary of some of the methods used to determine 
various atmospheric gas distributions is given below. Shown 
in Figure 1 are the experimentally determined mixing-ratio 
profiles (the mixing-ratio of a gas is the ratio of the 
number of gas molecules to the number of air molecules) of 
some atmospheric gases. The range of research being done 
on the atmosphere makes an exhaustive summary impossible. 
However, the following treatment illustrates the wide 
variety of experimental methods employed. 

In most in situ experiments stratospheric gas samples 
are collected at several altitudes by balloon, aircraft, or 
rocket. These samples are then analyzed by various methods. 
Mixing-ratio profiles of Cl, CIO, and OH have been deter- 
mined by laser induced fluorescence of samples obtained by 
balloon (4-6) and by aircraft (7,8) . Others have analyzed 
gas samples obtained from airplanes (9,10) and balloons ( Il- 
ls) by gas chromatography in order to determine the molecular 
abundances of CCI 2 F 2 , CC1 4 , CCl^F, CFCl^/ CF 2 CI 2 , CHCIF 2 , 

ch 3 cci 3 , ch 4 , co, c 2 ci 3 f 3 , c 2 ci 4 , c 2 hci 3 , c 2 h 2 , c 2 h 4 , c 2 h 6 , 

C 3 Hg, CgHg , C-yHg, H 2 / and ^0. Another method used to 
obtain abundances of CH 4 , CO, CO 2 , H 2 , *^0, ^0, and 0 3 was 
mass spectroscopy of samples obtained by a rocket (16) and 
a balloon (17) . 
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VOLUME MIXING RATIO 


Figure 1: The vertical distribution of several atmospheric 

gases after' Gille (1) . 
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Experiments which depend on the collection of gas 
samples are subject to the problems of contamination of the 
sample due to impurities on the container walls, contamina- 
tion due to leakage of the sample during storage, and the 
loss of active gases due to chemical reactions. Some 
investigators have avoided the problems associated with gas 
sampling by using dyes deposited by balloon or rocket which 
react chemiluminescently with the atmospheric constituent 
of interest. The spectral radiance of the sky due to these 
reactions is a measure of the concentration of the atmospheric 
constituent. This method has been used to determine the 
stratospheric profiles of NO (18) and 0^ (19-21) . 

In situ measurements of atmospheric composition allow 
only local regions of the atmosphere to be studied. A com- 
plete understanding of the atmospheric interactions requires 
global monitoring of the constituent distributions as well 
as the temporal variations of these distributions. Rapid 
observations with large spatial coverage have been made by 
remote-sensing techniques (22-38) . 

One remote-sensing technique which has been used is the 
collection and analysis of atmospheric emission spectra. 
Atmospheric emission spectra have been obtained from 
instruments at ground level (22) and in balloons (23). 

These spectra have been analyzed to obtain mixing-ratio 
profiles for CIO and OH. The radiation emitted by the * 
atmosphere has also been studied using satellite-borne 
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instruments and the "onion skin peeling" technique (24,25) 
to obtain the mixing-ratio profiles of CF 2 Cl 2 , H 2 0 ' HN0 3 / 
and NC> 2 . The radiative transfer problem for emission from 
atmospheric constituents is much more difficult to para- 
meterize than spectral absorptance. This limits the 
information obtainable from emission spectra analysis. 

Although adequate for monitoring minor constituents (ppm) 
and determining atmospheric temperature and pressure profiles, 
emission spectroscopy is generally not as useful as solar 
absorption spectroscopy for the study of trace gases (ppb) . 

Trace gases can be studied by solar absorption spectra, 
since most gases of interest have strong unobscured absorp- 
tion bands in the infrared, but accurate analysis of such 
spectra requires sophisticated analytical techniques and 
usually many simplifying assumptions. Many techniques have 
been developed for the interpretation of absorption spectra. 
Infrared absorption spectra obtained from the ground (26-28), 
balloons (29-32), and airplanes (33,34) have been used to 
determine the mixing-ratio profiles of CH^ , CO, COS, C0 2 , 

HCN, HC1 , HF, HN0 3 , H 2 CO, H 2 0, NO, N0 2 , N 2 0, and 0 3 . 

Visible absorption spectra obtained from balloons (35) have 
been used to determine the profile .of N0 2 . Ultraviolet 
spectra from balloons (36,37) and satellites (38) have been 
used to determine the stratospheric ozone profile. 

The global and temporal variations of many trace gases 
over a large altitude range can only be monitored by 
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absorption spectra obtained from satellites. However, few 
infrared absorption spectra have so far been obtained from 
satellite-borne instruments. The Atmospheric Trace Molecules 
Spectroscopy (ATMOS) experiment represents an attempt to 
fill this gap. ATMOS will allow the monitoring of many gases 
(some of which are listed in Table 1) up to 100 km above the 
earth's surface, nearly doubling the altitude range covered 
by balloon-borne instruments. The analysis described in 
this work has been done in support of the ATMOS project 
(however, the results obtained in this work can be used to 
aid the interpretation of many types of solar spectra) . 

A description of the ATMOS instrument and experiment are 
given below. 

We have attempted to determine the limits to the 
information retrievable from analysis of the spectral data 
which is expected to be obtained from this experiment as 
well as other experiments. Comparisons have also been made 
of the results expected from satellites to those obtainable 
from balloons and ground-based observers. 

1.2 THE ATMOS EXPERIMENT 

The ATMOS instrument, shown in Figure 2, contains 
a dual moving cat's-eye interferometer. The interferometer 
has been designed to make the instrument insensitive to 
angular and lateral motion of the moving components. This 
was accomplished by replacing the plane mirrors of a 
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Table 1: A compilation of some atmospheric gases expected 



to be detectable by the ATMOS 
1980) . 

experiment (as 

of 

Category 

_1 

Category 2 

Category 3 

CO 


HC1 

N 2 °5 

S0 2 

CN 

o 

u 


HF 

CH 3 C1 

H 2°2 

h 2 o 


NO 

nh 3 

ho 2 

n 2 0 


N0 2 

CH 3 F 

H0 2 N0 

ch 4 


CIO 

CH 3 Br 

CFC10 

HDO 


HN0 3 

F ll- 

cf 2 o 

°3 


OH 

F 12 

HOCl 




ciono 2 

no 3 

Category 

1 : Minor 

\ 

constituents 




Category 2: Trace constituents which have been observed. 

Category 3: Trace constituents of specific interest to 

ATMOS objectives which have not been quanti- 
tatively measured by remote sensing. 
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conventional Michelson interferometer with cat's-eye 
retroref lectors . In addition the radiation is double passed 
through the arms of the interferometer before being sent to 
the detector in order to eliminate distortions due to the 
lateral motion of the instrument. The HgCdTe detector is 
cooled by liquid nitrogen. 

The supporting equipment consists of a suntracker, 
a telescope, and a camera. The suntracker consists of two 
single-axis, motor-driven mirrors and a silicon diode 
assembly for sensing and controlling the position of the 
mirrors with respect to the sun. An f/2 telescope is used 
to focus the radiation received by the instrument into a 
beam suitable for the interferometer. A 16 mm' camera 
records the image of the sun on the entrance aperture to 
the instrument in order to determine the portion of the disc 
viewed, to insure, if necessary, that data obtained for 
a field of view containing a sunspot may be identified and 
excluded from the data analysis. 

The ATMOS instrument is scheduled to fly on Spacelab 3 
in 1984. The experiment to be conducted involves viewing 
the sun with the ATMOS instrument at periods when portions 
of the earth's atmosphere occult the sun. Measurements will 
be made of the incident solar radiation as a function of the 
optical path difference of the beams in the instrument. 
Fourier transforms of these collected interferograms will 
give solar absorption spectra covering the region from about 
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2 pm to 16 pm. Measurements obtained for various radiation 
paths through the atmosphere will allow the determination of 
the mixing-ratio profiles of atmospheric constituents. 
Because of the orbital geometry of the space shuttle and the 
rate at which the interferograms are measured by the ATMOS 
instrument as well as the broad wavelength coverage, 
simultaneous mixing-ratio profiles with good vertical reso- 
lution of many constituents will be obtained together with 
the global variations of these profiles. This will give 
a three-dimensional view of the constituent abundances, 
a view which is crucial to an understanding of atmospheric 
processes. Some of the molecules expected to be detectable 
by ATMOS are listed in Table 1. Little is known about the 
global distribution of many of these gases, many of which 
are important in the reduction of maintenance of the ozone 
layer. 

1.3 SUMMARY OF THIS WORK 

In this work an attempt is made to determine the amount 
of information about the distribution of absorbing gases in 
the atmosphere obtainable from analysis of the data to be 
obtained from the ATMOS experiment. The results obtainable 
for satellite remote-sensing have also been compared to the 
more common remote-sensing balloon-borne and ground-based 
experiments . 
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One method which has been used by others to analyze 
spectra obtained by satellite-borne instruments is the 
"onion skin peeling" technique (24,25). In this technique 
the radiative transfer problem is modeled in the many 
atmospheric layers traversed by the radiation one at a time. 
Thus errors introduced by the analysis method into the 
solution in previous layers will be compensated for in 
successive layers, propagating the errors throughout the 
obtained mixing-ratio profile. In this work we have 
developed a technique which analyzes all portions of the 
radiation paths of several solar rays simultaneously assuming 
no a priori knowledge of the shape of the profile to be 
retrieved. Thus any error in the retrieved mixing-ratio 
profiles is due only to the experimental noise and uncertain- 
ties in the experimental conditions and not to inadequacies 
in the analysis technique. The results in this work repre- 
sent the optimum obtainable for the specified experimental 
conditions . 

In Chapter 2 the uncertainty and the vertical resolution 
as a function of altitude obtained in the retrieval of 
a vertically-co.nstant , mixing-ratio profile are estimated by 
analyzing synthetic, equivalent-width data using a linear 
least-squares technique. Estimates of the uncertainty and 
vertical resolution are obtained for observing instruments 
at satellite and balloon altitudes as well as for ground- 
level observations. 
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In the third chapter the effects of systematic errors 
in the parameters describing the instruments orientation 
with respect to the earth as well as those describing the 
thermal structure of the atmosphere on the retrieval of 
mixing-ratio profiles are investigated using a linear 
least-squares technique. A method is described by which 
many of these systematic errors can be eliminated in the 
analysis . 

In the fourth chapter the information contained in the 
shape of spectral lines is estimated. Spectral data obtained 
from a recent ground-level test of the ATMOS instrument are 
analyzed using a non-linear least-squares technique to 
retrieve information about the mixing-ratio of N 2 0. 

in this work we have not attempted to model the 
smearing effects on spectral data due to the motion of the 
observing instrument (39) or the finite field of view of the 
instrument (40). These effects (and others) need to be 
considered in order to accurately determine the information 
content of such data. Although the modeling of such smearing 
effects will modify the results obtained in this work, we 
have presented a method of analysis which allows the deter- 
mination of this information content as well as the structure 
of retrieved mixing-ratio profiles. 
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Chapter II 

THE RESOLUTION OF RETRIEVED MIXING-RATIO PROFILES 
2.1 INTRODUCTION 

The analysis of solar spectral data to retrieve mixing- 
ratio profiles of atmospheric constituents requires assump- 
tions to be made about the kind of information retrievable 
and the quality of this information. 

If it is assumed that a continuous mixing-ratio function 
of known shape can be retrieved and the shape of the assumed 
function is incorrect or if the data contain insufficient 
information, then it is possible to obtain a profile with no 
apparent systematic errors but very large uncertainty. 

The inclusion of the a priori information (1) forces a 
solution which may not be correct. If it is assumed that 
only independent, discrete values of the mixing-ratio 
profile can be obtained, it is necessary to know how many 
values are obtainable from the data, their spacing, and their 
uncertainty in order to determine the appropriate analytical 
method. 

The extraction of geophysical information from gross 
earth data containing noise has been considered by Backus 
and Gilbert (2-7) , A method was derived in which the 
vertical resolution of retrieved geophysical profiles was 
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found in terms of the uncertainty in the retrieved profile. 
The method of Backus and Gilbert' has been used to estimate 
such varied quantities as the vertical resolution of 
atmospheric temperature profiles (1,8,9), aerosol size dis- 
tributions (10), and the conductivity in the earth's mantle 
(11) . No such studies have been done on the vertical 
resolution of constituent mixing-ratio profiles obtained 
from absorption spectra. 

In this chapter the problem of obtaining a mixing-ratio 
profile from spectral data containing noise by assuming all 
the other physical parameters describing the experiment are 
exactly known is investigated. The data analyzed consisted 
of absorber amounts, calculated by using Eqs . (A. 95) - (A. 97) 

(See Appendix A) , for various paths through the atmosphere 
as shown in Figure 3. These absorber amounts, as shown in 
section 2, could be obtained from equivalent width analysis 
of weak lines appearing in the spectra. 

In section 3 the method of Backus and Gilbert is applied 
to the extraction of mixing-ratio profiles from absorber 
amount data. The analysis shows that only discrete mixing- 
ratio values are retrievable and the uncertainties in the 
retrieved values are minimized if the sum of the square 
residuals (RSS) , where the residual is the difference 
between the measured absorber amount and that obtained using 
the retrieved profile, is a minimum. It was found that the 
uncertainty increases and the resolution degrades with 
altitude . 
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Figure 3: The path of radiation through a layered atmosphere 

where each layer is characterized by lower geo- 
potential altitude Hj_, upper altitude and 

temperature Tj_ at altitude Hj_ . The optical 
airmass encountered by the ray in element of path 
ds is da. 
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' An analytical method which uses the minimization of the 
RSS (described in section 4) was developed to retrieve 
mixing-ratio profiles from absorber amount data. The 
atmosphere , as shown in Figure 3 , was divided into concen- 
tric shells at which the temperature, pressure, and density 
are defined. Absorber amounts were calculated for each ray 
path by assuming a constant mixing-ratio profile for the 
absorber. Random noise was added to the absorber amount 
data sets. The mixing-ratio values were then retrieved from 
these simulated data sets by using P9R (12) , a linear least- 
squares computer program. 

This program retrieved a set of local mixing-ratio 
values and their standard deviations. The standard devia- 
tions were used to estimate the uncertainty in the retrieved 
profiles. In section 5 a general relation is given for this 
uncertainty in terms of the altitude and the mixing-ratio 
value spacing. 

The vertical. resolution was determined in the sixth 
section by using the uncertainty function obtained in the 
previous section. The forms of these expressions agree with 
those predicted by the Backus and Gilbert method. 

The mixing-ratio value spacing needed for a particular 
analysis is determined by the vertical resolution. This 
spacing should not be less than half the vertical resolution 
at any altitude and thus it also determines the number of 
retrieval values. By reducing this number the amount of 
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computer time required is reduced, an important factor when 
the retrieval is accomplished by modeling the spectral line 
shapes of several overlapping lines. In the seventh section 
the appropriate retrieval parameter spacing for two values 
of the SNR is discussed. 

Finally it is desirable that the derived relations for 
the uncertainty and the vertical resolution be extended to 
mixing-ratio profiles of arbitrary shape. In the final section 
comparisons were made of retrievals from data obtained from 
other profiles. The functional forms obtained for the 
uncertainty and the vertical resolution in the fifth and sixth 
sections were found to be adequate to describe that obtained 
for these profiles. 

In brief the goal of this chapter is to determine the 
maximum amount of information about the mixing-ratio profile 
retrievable from absorber amount data and the quality of the 
retrieved values, with the object of constructing the optimum 
retrieval scheme for a particular data set. 

As an aid to the reader, the symbols used in this work 
are listed in Appendix B. 

2.2 EQUIVALENT WIDTH ANALYSIS OF SPECTRA 

In this chapter we investigate the analysis of absorber 
amount data. A set of solar spectra obtained from solar 
radiation passing through various portions of the atmosphere 
contains information about the distribution of absorber along 
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the path traversed by the radiation. The total absorber 
amount along the ray can be determined by equivalent width 
analysis of the spectra. 

Let I v be the intensity of radiation per unit area from 
a nearly parallel beam. After passing through a path- 

length da m of absorber, the change in intensity dl v is 

dl v = -g v k v I v da m (la) 

where g v is some function representing the effect of the 
instrument on the observed radiation and k v the absorption 
coefficient at frequency v. The total change in the energy 
in the beam dl is 


dl 


oo 



Ida m dv 


(lb) 


If the instrument function is defined only over some fre- 
quency interval Av the total change in energy is 

■- / 


dl 


Av 


gykvlda^dv 


(lc) 


and, if Av is small* (i.e., high spectral resolution) then 
k v and I are constant over the integral or 


dl = -k v I v da m J g v dv 
Av 


/- 


(Id) 


Thus the total observed change in energy is scaled' by the 
integral of the instrument function. Normalization of 


*Nielsen et al . (13) have shown that the equivalent 

width is independent of the instrument function width. 
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Eq. (Id) gives 


dl = -k v Ida m 


(le) 


or upon integration 


I (v) 


= ex P [ ”k v a m ] . 


(If) 


If the transmittance T v is defined as I(v)/i Q , then 


T v = exp[-k v a m ] . 


(2a) 


The spectral absorptance at frequency v is defined by 


\ = 1 - T v • 


(2b) 


For a single spectral line the area under the absorption 
curve shown in Figure 4 is given by 


W(a m ) 


■/: 


A,,dv . 


( 3 ) 


The area W(a) is equal to the area under a rectangle of 
width W(a) as shown in Figure 5 and W(a) is defined as the 
equivalent width of the spectral line of Figure 4. Combining 
Eqs. (2 ) and (3) gives 


W(a m ) 


-L 


(1 - exp[-k v a m ] )dv . 


(4) 
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Figure 5: The equivalent width W of a single absorption 

line. The area of the rectangle is equal to the 
area between the absorption line and T =1. 
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If k v a m <<l then the exponential in Eq . (4) can be 

approximated by the first two terms in the series expansion. 
Thus the equivalent width becomes 


W(a m ) 



k v amdv . 


(5) 


Since ajjj does not depend on the frequency. 


w ( a m) = a m / k v dv 


( 6 ) 


The intensity S of the spectral line is defined as 


■I 


S =/. CO k V dV ' 


(7) 


so 


W(a m ) = Sa m . (8) 

For weak lines (A V <<1) the equivalent width is propor- 
tional to the absorber amount along the path provided S is 
not dependent on the temperature. If the line intensity is 
known and the area under the absorptance curve is measured, 
then the absorber amount if found by 

am = W(a m )/S . (9) 
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This expression was obtained for a homogeneous gas 
sample. The atmosphere is inhomogeneous, and the temperature 
and density vary along the path. The absorption coefficient 
k v is temperature and pressure dependent and the intensity 
S is temperature dependent. The problem of the inhomogene- 
ous path is complicated. However, if the atmosphere is 
divided into concentric layers the ray path can be divided 
into a series of atmospheric portions which can be considered 
to be homogeneous. Provided A V <<1, since the equivalent 
widths of weak lines are additive, the total absorber amount 
can be estimated for atmospheric paths. 

For larger values of A v Eq. (9) is no longer true. In 
this case the expressions for W(a m ) are more complicated 
(14,15) and it is more difficult to extract information about 
the absorber amount. We have considered only the case of 
weak spectral lines in this chapter. This represents not 
only the simplest set of data to analyze but also, because of 
the additive nature of the equivalent widths, it yields the 
maximum information. 

The equivalent width method is useful when the true 
spectral features are narrower than the spectral resolution, 
since the equivalent widths obtained from actual spectra are 
essentially independent of the shape of the instrumental 
response function and its width (15). Equivalent width 
analysis has been used by many groups including Buijs et al 
(1.6) who determined the mixing-ratio profiles of HC1 and HF 
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between 15 and 30 km altitude from spectra obtained by 
a balloon-borne instrument. Other investigators (17,18) 
have used equivalent width analysis to obtain mixing-ratio 
profiles. for CH 4 (17) and H 2 CO (18) from ground-level 
observations. In the two ground-level retrievals a priori 
assumptions were made about the shapes of the mixing-ratio 
profiles to be obtained. Even so the vertical resolution 
was very poor. Some of the reasons for this are discussed 
below. 

2.3 THE THEORY OF BACKUS AND GILBERT 
2.3.1. Genera l formalism 

The total constituent absorber amount a m encountered by 
a solar ray traversing the atmosphere is given by 

m(r)da ( 10 ) 

where m(r) is the constituent mixing-ratio by volume of 
absorber m, da is the number of atmospheric molecules per 
unit area in path element ds and r is the distance from 
earth's center (see Figure 3). The retrieval of the mixing- 
ratio profile from a set of absorber amounts involves the 
inversion of Eq. (10). The inversion of integral equations 
of this type is often required in atmospheric problems and 
these inversions are often difficult. If no a priori assump- 
tions, are made about the functional form of m(r), only 
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discrete, mean values of m(r) can be obtained at values of 
r corresponding to the peak in da for each absorber amount. 
Also the noise in the experimental absorber amounts causes 
uncertainties in the discrete mixing-ratio values retrieved. 

Backus and Gilbert (2-7) have considered the inversion 
of equations of this type and have shown the noise causes 
an unavoidable trade-off between the vertical resolution and 
the uncertainties in the derived profiles. 

Consider N 0 measurements of the absorber amount along 
N 0 different ray paths through the atmosphere. The mixing- 
ratio profile m(r) is obtained by simultaneous solution of 
the N 0 linear, first-order, integral equations 

/■ r TOP dai 

a mi ( m ) = / m(-r) dr (11) 

^ r Ti 

where r T ^ is the radius of the sperical shell which contains 
the closest approach or tangent height of the ray and r TOp 
is the radius of the spherical boundary of the top of the 
atmosphere . 

It is convenient to introduce the dimensionless para- 
meter 


r - r 


Y = 


Ti 


r TOP " r Ti 


In terms of this parameter Eq. (11) becomes 


( 12 ) 
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-gpr dy 


Thus the N 0 a m; j_ (m) measurements are given in terms of the 
dai(y)/dy weighting functions. 

A composite weighting function can be constructed from 
the N Q weighting functions of Eg. (13) by defining 


D (y) = Z b ± 
i = l 


dai (y) 


where the constants b]_,...,b N are chosen such that D(y) is 

1N o 

unimodular or 


D(y)dy = 1 . 


We would also like to choose b]_,...,b No so that D(y) is 
weighted heavily near a value of y Q where the local average 
<m e >y 0 of m e (y) is desired. The function m e (y) is the true 
mixing-ratio profile to be retrieved. It is desirable to 
have weighting function D(y) closely approximating the 
delta function about y Q . Thus the average <m e > which is 
retrieved is nearly equal to the true value m e (y D ) . 

The spread of weighting function D(y) from y was defined 
by Backus and Gilbert (7) as 


S(y 0 , D) =12/ (y - y Q ) 2 D(y) 2 dy 
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This spread is a measure of the width of the peak in the 
weighting function D(y). This width approximates the pos- 
sible vertical resolution in the retrieved profile. 

An inner product in the function space of the kernels 
of the integral equations given from Eq. (13) can be defined 


<m, D > 



m(y)D(y)dy . 


(17) 


The variation of the retrieved profile from the true mixing- 
ratio profile m e (y) is defined as 


A<m e , D> = <m, D> - <m e , D> 


(18) 


In terms of the inner product 


A<m e , D> 


N o 

Z b: 
i=l 


/ 


1 

.[m (y) - m e (y) ] 
o 


da^ (y) 
dy 


dy . 


(19) 


The difference between measurement i of the absorber amount 
and the absorber amount given by the true mixing-ratio pro- 
file m e (y) is 


AE 




[m (y) - m e (y) 3 


daj_ (y) 
dy 


dy 


( 20 ) 


Thus the variation is 


A<m e , D> 


E b^AEi . 
i=l 


( 21 ) 
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0 


Since AEj_ is the noise in measurement i, then AEj_ = 
if the noise is randomly distributed over the N 0 data values. 
The square of the variation is 


(A<m e , D> ) 2 


N 0 N q 

( Z b i AE i ) ( E b • AE . ) 
i=l j=l 3 3 


( 22 ) 


The mean square variation is given by 


(A<m^, D > ) 


N, 


I b i bjAE i AEj 
i, j=l 


(23) 


In terms of the matrix formulation 


(A<m p , D>) 2 .= b T Eb , 


(24) 


where E is the covariance matrix of the measured values a m i . 
The ith mean square variation is equal to the standard 
deviation a-j_ of the ith retrieved value (16) so 


2 U T_. 
a = b Eb . 


(25) 


The optimum combination of retrieval value uncertainty 
and vertical resolution is found by minimization of the 
linear combination 


Q (y) = qS (y) + (l-q)acr 2 (y) 


(26) 
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2 

where a is a factor which converts a into the same units 
as S (y) . The relationship obtained for S(y) in terms of 
a(y) is called the trade-off curve. The appropriate combi- 
nation of spread and standard deviation can be found for 
all values of 0 <_ q < 1. 

2.3.2 The theory applied to mixing-ratio retrieval 

We wish to determine the best possible uncertainty 
(standard deviation) and vertical resolution (spread) in 
the retrieval of mixing-ratio profiles from absorber amount 
data. The smallest uncertainty corresponds to q = 0 and the 
highest vertical resolution to q = l in Eq. (26). These 
represent bounds for a and S. 

The kernel da^/dy of Eq. (14) is proportional to the air 
density (y) along the ith path through the atmosphere. 

Since the composite weighting function is required to be 
unimodular the coefficient of the ith kernel is given by 

b i = 1 . (27) 

f 1 da i 

Jo sr d Y 

For an isothermal atmosphere 

p (y) = p (0) exp[-y/H s ] , (28) 

where H s is the scale height given by 
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( 29 ) 


R*T . 
s Mg ’ 


Thus 

bjLaexp[y/H s ] . (30) 

If the N 0 measurements of are independent the cross 
terms of the covariance matrix E are zero, thus 

■ E = ct e 2 T , (31) 

where a E is proportional to the noise in the experimental 

2 

data. From Eq. (25) a becomes 
2 2 T 

a = a E b l b (32) 


or 



(33) 


Finally, combining Eqs. (30) and (33) we get 


o a 


Y_ 

e H s 

SNR 


(34) 


since is proportional to the experimental noise. Using 
the definition for y and the relation 
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r 


R e + A 


(35) 


where A is the geometric altitude above the earth's surface, 
Eq . (34) becomes 


a a 


A 

H_ 


SNR 


(36) 


Thus the minimum uncertainty in the retrieval increases 

exponentially with altitude. 

The other limit that is needed is the minimum vertical 

> 

resolution (spread) . Rewriting Eq. (17) in terms of the 
composite weighting function gives 

N o ..1 2y 

Sa E b i b i / (y-y 0 ) 2 e H s dy . (37) 

i,j=l J J o 

Upon integration Eq. (37) reduces to 

«o -f- 

Sa E b i b i e Hs (38) 

ifj=l 


or 



i/ j=l 


b i b j 



Thus the spread reduces to 


Sa 


H- 
e s 


(39) 


(40) 
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or, in terms of the altitude 




Sa- 


(4i; 


Therefore it is expected that the minimum uncertainty 
will increase and the optimum vertical resolution will 
degrade exponentially with altitude. The constants of pro- 
portionality in Eqs . (36) and (41) cannot be determined by 

the Backus-Gilbert method alone. They are obtainable from 
simulated retrievals on synthetic absorber amount data. 

This theory does allow the method of analysis which will 
give a minimum uncertainty to be determined. In the following 
derivation the proper method of analysis is determined. 

2.3.3 Criterion for minimum uncertainty 

The integral relations in Appendix A, Eqs. (A. 95) - (A. 97), 
for the optical airmass encountered by a refracted ray can be 
approximated by numerical integrations if the atmosphere is 
divided into concentric layers in which the relations (A. 95)- 
(A.97) can be applied successively. In this formulation the 
jth absorber amount value is given by 

N o 

a mj = ^ m i a ij ' (42) 

i=l 


where mj_ is the local mixing-ratio value (assumed constant) 
in layer i and a^j is the optical airmass encountered by the 
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jth ray in layer i. 

If it is assumed that the atmosphere is uniform in all 
variables except altitude (i.e. , the atmosphere is "horizon- 
tally" uniform) , then the N 0 weighting functions in Eq. (14) 
are identical and the composite weighting function becomes 


D(y, - b . 


(43) 


Using Eq. (43) the standard deviation given by Eq. (24) 
reduces to 


, N o 

fc> £ AE ± AE j 

i,j=l 


(44) 


Expansion of the j summation gives 


N r 


2 o w 2 

o = b^ E AE-; + AE-, AE 0 + ...+ AE m AE. t . 

i=1 1 12 N o Nq -1 


(45! 


If N 0 is large then 


N r 


AE X (AE 2 + AE 3 + . . . + AE n ) lAEi . 2 AE i ) , 


i=l 


(46) 


so Eq. (45) becomes 


Nr 


N, 


a 2 = b 2 E (AEi) 2 + AE x E AE i + 
i=l i=l 


N o 

. + AE^ E AEj_ 
i=l 


(47) 
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But 


No 

Z AE]_ - 0 , (48) 

i=l 

if the noise is randomly distributed over the N 0 data values 
Thus the mean square deviation is given by 

N 

a 2 = b 2 Z ( AE j_ ) 2 (49) 

i=l 

or, in terms of the measured absorber amounts, a m . 

N 0 

•o' 2 - b 2 * . (50) 

i=l e 

The residual sum of squares (RSS) is defined by 
N 

° 2 

RSS - E ( a mi ~ a m e i^ • (51) 

i=l e 

Thus the uncertainty a is minimized if, and only if, the RSS 
is minimized. In order to analyze synthetic absorber amount 
data we chose a computer program (P9R) which performs 
linear inversion by RSS minimization. This guaranteed that 
the uncertainties in the retrieved values would be minimized 
The sections which follow describe the method used to 
calculate the synthetic absorber amounts and the method of 
retrieval of the mixing-ratio profile. These have been used 
to find the proportionality constants in Eqs . (36) and (41) 
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from analysis of synthetic absorber amount data calculated 
assuming a constant mixing-ratio profile. These constants, 
in general, will depend on the earth-observer geometry as 
well as the mixing-ratio profile retrieved. It is useful to 
determine these constants so that the magnitudes of the 
uncertainty and the vertical resolution can be estimated. 

The results of this analysis are given in sections 7 and 8. 

2.4 RETRIEVAL METHOD 
2.4.1 General comments 

We have shown that the total absorber amount along 
a ray path could be obtained from the equivalent widths of 
a single spectral line in spectra corresponding to different 
paths and that the uncertainty in the mixing-ratio profile 
retrieved from such data is minimized if the RSS is minimized. 
If the retrieval problem is represented by Eq. (42) , the 
mixing-ratio values can be retrieved by linear inversion. 

The path of radiation through the atmosphere can be 
characterized by the path length segments and the direction 
of travel through concentric layers or by the optical air- 
mass encountered by the ray in each layer. The airmass can 
be obtained from Eqs . (A. 95) and (A. 97) if the temperature 

and air density at each layer boundary are known. If the 
airmass elements aj_j of Eq . (42) are calculated for each 
segment of the N Q ray paths, Eq. (42) can be inverted to 
give up to N Q local mixing-ratio values. 
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In the following section the method used to construct 
an atmospheric model, necessary for the application of 
Treve ' s formalism (19,20), is discussed. Next the ray path 
calculation technique (i.e., the calculation of the a^j 
values), using Eqs . (A. 9 5) - (A. 97) of Appendix A, is given. 

Finally, the linear least squares method of retrieving 
local mixing-ratio values from a set of total absorber amount 
values is presented. 

A more detailed presentation of the atmospheric model 
and the optical airmass calculation is given in Appendix A. 

It may be helpful for the reader to examine Appendix A 
before continuing. 

2.4.2 Construction of the atmospheric model 

The optical airmass encountered by radiation passing 
through the atmosphere was obtained using the relations 
derived in Appendix A. The complicated nature of the 
integrands required that the integrations be performed 
numerically. In order to do this the atmosphere was 
divided into concentric spherical layers, as shown in 
Figure 3, characterized by the heights of the upper and 
lower boundaries or break points of each layer and the 
temperature and density at each break point. If the height 
and temperature of each break point, together with the 
pressure at the bottom of the lowest layer are given, the 
value of each of the quantities in Eqs. (A. 95) - (A. 97) can be 
determined. 
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If the latitude of observation is known the geopotential 
height of each break point can be obtained from Eq . (A. 14) . 

The pressure at each break point can be calculated from 


9o m 


R*L 


i+l 


= 


1 T i + L i (H i+1 -H i ) 


g 0 M 

^i -rSt" ( H i+1- H i } 


1 Li ^ 0 

if 


L i = ° 


(52) 


where Li, the temperature gradient in the ith layer, is 
given by 


T 

i+l 


H 


i+l 



C53) 


The mass density at each break point was found using the 
ideal gas law 


Mp i 

p mi = r*tT 


(54) 


The number density was obtained from 
N A P i 

Pvi ~ R*Ti ' (55) 

The index of refraction of air for T = 288 K and P = 1 atm 
was given by Edlen's formula (21) 


(n-l)xlO 8 = 6 43 2.8 + 2949 810.(146 - 1/X 2 )" 1 
+ 25540. (41 - l/l 2 ) _1 , 


(56) 
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where ^ (ym) is the wavelength of the observed radiation. 
The index of refraction at each break point was found from 
the Lorenz-»Lorentz formula (22) 


2 -i 

i n i+i - 1 
P vi+1 n? +1 + 2 


p vi n? + 2 


(57) 


Thus far we have constructed an atmosphere consisting of 
layers defined by the various physical parameters at the 
boundaries. Within each layer the pressure, temperature, 
and index of refraction varied according to Eqs . (52), (55), 

and (57) respectively. The temperature gradient L-^ was 
assumed constant throughout a given layer. 

We assumed that the atmosphere was "horizontally" uni- 
form. The temperature and density at altitude A were assumed 
independent of the position defined by the latitude and 
longitude of the earth. The size of the errors which result 
from this assumption are considered in the next chapter. 


2.4.3 Calculation of the optical airmass 

Now that the atmospheric model is defined the airmass 
encountered by a ray traversing a refracted path through the 
atmosphere can be computed. 

Consider radiation with wavelength A reaching an instru- 
ment at height A 0 at an apparent zenith angle Z Q . If we 
assume that the source of the radiation is an infinite 
distance away from the earth and that A, A Q , and Z G are 
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known, the path of the radiation and the optical airmass 
in each layer are obtained from Eqs . (A. 95 ) - (A. 97 ) . 

The optical airmass values were computed using a version of 
the program described by Snider (23) and Snider and Goldman 
(24) modified to allow observer heights up to 500 km. The 
method used for the numerical integration was Gauss quadra- 
ture with K = 6 (25) . Treve (20) has shown that this approxi- 
mation is accurate over the lowest 90 km of the atmosphere. 
Since all the atmosphere of interest in our investigation was 
below 60 km this was considered adequate for our purposes. 

There are two classes of rays through the atmosphere. 

The first reaches the observer at some apparent zenith angle 
Z o £90°. This case is the simpler of the two, since the 
observer height is the lowest penetration of the ray in the 
atmosphere. The second case is Z o >90°. For these apparent 
zenith angles the ray passes through a portion of the atmos- 
phere below the observer. However, if the closest approach 
to the earth of the ray (the tangent height) can be deter- 
mined, the ray path can be divided about the tangent height 
into two paths of the first type. 

2.4.4 The two possible geometries 

If Z o <_90°, is shown in Figure 6, the angle of entry 
Zi in layer i is related to the angle of exist Z Q j_ by 

Zi = Wi + Z Q i - «i (58) 
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(see Figure 77 in Appendix A) where the angle is found 
from Eq. (A. 96) and from Eq . (A. 95). The angle of entry 

Zj_ is equal to the angle of exit from the layer above. 

Thus the path of the ray and the total airmass encountered 
can be obtained by successive used of Eqs . (A. 95) - (A. 97) from 

the layer of the observer to the top of the atmosphere (of 
the source height A s ) . The connection between each layer is 
provided by 

z i = z oi+l * (59) 

The more difficult case of Z o >90 ( ? is shown in Figure 7. 
In this case the ray traverses a portion of the atmosphere 
below the observer height A 0 . The closest approach of the 
‘ray to the earth is called the tangent height A T . This ray 
path can be divided into two portions; that from the tangent 
height to the source and that from the tangent height to the 
observer. Since the apparent zenith angle at the tangent 
height is Zt = 90° the problem reduces to two problems of the 
first type discussed. 

The tangent height can be determined by use of Fermat’s 
principle (see Eq. (A. 45)) 

nrsin Z = n 0 r 0 sin Z Q . (60) 
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Figure 7: The slant path geometry for a satellite observer 

and a ray reaching the observer at apparent 
zenith angle Z o >90°. The closest approach of the 
ray to the earth (the tangent height) is A^ . The 
ray path is divided into two portions CB and BA 
about the tangent height. 
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At the tangent height sin(Z T ) =1, so Eq. (60) reduces to 

nr = n 0 r 0 sin Z Q . (61) 

The tangent height is found by solving the transcendental 
equation 

nr - n 0 r 0 sin Z Q = 0 , (62) 

where r is obtained from Eq . (A. 14) 



Y 


and H t is the geopotential equivalent of A T , and y is the 
conversion factor from geometric to geopotential meters. 

As noted in Appendix A all altitudes used in the airmass 
computation are in geopotential meters. At the earth's sur- 
face the difference between geometric and geopotential 
altitude is 0 km. At 50 km the difference is about 0.4 km. 
Thus, if geometric altitudes are used instead of the geopo- 
tential altitude the error in tangent height below 50 km can 
be as large as 0.4 km. 

Once the tangent height is determined the ray in 
Figure 7 can be replaced by the two ray segments BA and BC 
and the problem reduces to solution of the ray tracing 
problem from the tangent height to the source and from the 
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tangent height to the observer. The total airmass is found 
by successive application of Eqs . (A. 95) - (A. 97 ) in each of 

the layers in the two portions of the ray. 

Thus the airmass elements a-j_j in Eq. (14) can be calcu- 
lated for any set of ray paths. Synthetic absorber amount 
data can be obtained from Eq. (14) if a mixing-ratio profile 
is assumed. Now that the a^j matrix and the M 0 a m values are 
known Eq. (14) can be inverted to give the local mixing-ratio 
values. The retrieval method used in this work is described 
below. 

2.4.5 Mixing-ratio profile retrieval method 

The synthetic amount of absorber along a ray path is 
given in terms of the local mixing-ratio values by 

N 

a m (Z oi ) = l n^a^ + e ± (64) 

j = l 

where Z 0 j_ is the apparent zenith angle of the ith ray and 
is the noise contained in the ith absorber amount. 

The formulation of Eq. (64) has implicitly assumed that 
there is no "horizontal" variation in the vertical mixing- 
ratio profile along the ray path. If, as in Figure 8, the 
mixing-ratio profile in the first half of the path is dif- 
ferent from that in the second half then Eq. (64) becomes 

N N 

a m( z oi ) = ^.^lj^j +Js j £ 1 m 2j a ij +e i ' (65) 
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DARK SIDE SUN SIDE 


Figure 8 : The case of an absorber with a mixing-ratio profile 

which is different for the two halves of the ray 
path . 
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which reduces to 


N 

a m (z oi ) = r + m 2 j)a i j + e ± . (66) 

j = l 

In this case, as shown in Figure 9, the mixing-ratio profile 
retrieved is the average of the two profiles. Such varia- 
tions in the mixing-ratio profile may occur for molecules 
such as NO or CIO (26,27) whose concentrations change 
rapidly at sunrise or sunset. In these analyses we have only 
considered the case of a constituent mixing-ratio profile 
which is "horizontally" constant. 

By making the simplifying assumptions discussed above 
the local mixing-ratio values in Eq. (64) can be retrieved 
with linear regression methods. The method used in this 
work (computer program P9R (12)) minimizes the RSS . As 
shown in section 3 this linear least squares technique 
retrieves local mixing-ratio values with minimum uncertainty. 
This program also has the capability to determine the number 
of values retrievable from a data set containing insufficient 
information. This feature was particularly useful for data 
obtained for a ground-level observer. 

The data sets analyzed by this method consisted of 125 
absorber amounts corresponding to 125 different apparent 
zenith angles Z Q j_ . If these are known the a^j elements in 
each layer can be obtained. For the purpose of the airmass 
computation we divided the atmosphere between 0 and 50 km 
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9: A mean "retrieved" profile for an absorber with 

a mixing-ratio profile which is different for the 
two halves of the ray path. 
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into 200 layers, each 0.25 km-thick. As shown in the next 
chapter this thickness is sufficiently small to avoid 
significant systematic errors in the airmass calculation. 

If only 125 absorber amount values exist in each data 
set, at most 125 local mixing-ratio values can be retrieved. 
The mixing-ratios in the atmosphere between the retrieved 
values were assumed to vary linearly with height. Thus if 
the local mixing-ratio value m^ in layer i is to be retrieved 
(see Figure 10) , the local 'mixing-ratio value in layer i+1 
is given by 

( m k+l ~ m k) 

mi+l = m^ + , (67) 

s 

where s is the number of atmospheric layers between the local 
mixing-ratio values m^ and m^+g to be retrieved. If the Nr 
mixing-ratio values to be retrieved are equally spaced with 
altitude Eq. (64) can be written as 

Nr s (2. - 1) (m k+ g - mj^) 

a m( z oi) = ^ ^ m k + i a i- , t (k - 1) s + &] 

k=l A=1 

+ eg . (68) 

This equation can be rearranged to give a linear combination 
of the local mixing-ratio values. This linear combination is 


51 




Figure 10: The structure of the retrieval atmosphere. 

Discrete mixing-ratio profiles were retrieved 
for layers k and k + 1. The ray path was traced 
through the S layers between k and k + 1. 
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a m(Z 0 i) - m]_ Z 
£=1 


S s + 1 - £ ^ r s £ 

~ a, £ + E m k E ^ a. r iBxll 

s 1 k=2 kj.^ s i, [ (k - 2) s + £ + 1] 


, S " 1 s - £ 

+ E a . 


£=1 


s i , [ (k - 1) s + £ + 1] 


] +E i • 


(69) 


Equation (69) can be written as 


Nr 

a m( z oi) = E m k a ik + e i ' 
k=l 


(70) 


where 


l il 


“ s + 1 - £ 

u — — — — d, , • 

£=1 S l£ 


(71) 


and 


^ rt A/ 

a ik = “S’ a i,[(k-2)s + £ + l] 

s-1 

_ s - £ 

+ Z a . 


£=1 


s i , [ (k - 1) s + £ + 1 ] , for k^l 


(72) 


Thus for any number Nr (<N q ) of local mixing-ratio values 
the inversion reduces to that of Eq. (64) . 

Since the N R equally-spaced, local, mixing-ratio values 
were assumed to be independent no a priori assumptions about 
the shape of the profile to be retrieved were necessary. The 
retrieval method (P9R) obtained both the local mixing-ratio 
value 1 % and’ its standard deviation . 
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This retrieval method was used to analyze absorber 
amount data obtained assuming various mixing-ratio profiles 
and observer altitudes. The results of these investigations 
are presented in the remaining sections of this chapter. 

2.5 THE UNCERTAINTY IN A PROFILE RETRIEVAL 
2.5.1 General comments 

We have shown in section 3 that the uncertainty in a set 
of retrieved mixing-ratio values is minimized if the RSS is • 
minimized. In section 4 a retrieval method utlizing RSS 
minimization was described. In this section this retrieval 
method is used on synthetic absorber amount data obtained 
assuming a constant mixing-ratio profile, shown in Figure 11, 
in order to estimate the minimum uncertainty in the retrieved 
mixing-ratio values. 

We have considered three typical observer geometries 
which represent the three types of occultation experiments. 
The first was the satellite geometry shown in Figure 12 A) 
where the observer height was 500 km. The second typical 
geometry is the balloon geometry of Figure 12 B) where A was 
30 km. The last geometry considered was the ground-level 
observer of A = 0 km. For the satellite observer all the 
Z Q i values were greater than 90°. For the ground observer 
all the Z Q values were less than 90°. The balloon observer 
is a combination of the satellite and ground geometries. 
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Figure 11: A mixing-ratio profile constant with altitude. 
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Figure 12: Occultation geometries for experiments performed 

at satellite A) and balloon B) altitudes. The 
observer position is given by 0. 
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Table 2: The distribution of the apparent zenith angles 

corresponding to the 125 data values computed for 
the 0 km observer. 


Angular 

range 

Angular spacing 
of adjacent rays 

Number of rays 
in each range 

10° - 

CTi 

O 

o 

2.0° 

25 

G\ 

O 

0 

1 

00 

Ul 

o 

0.5° 

50 

1 

o 

LD 

00 

O 

O 

0.1° 

50 
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Table 3 : Standard atmospheric values used for the retrieval 

of mixing-ratio profiles. 


Parame ter 

Assumed value 

latitude (which gives 
R e and g Q ) 

o 

in 

observer altitudes 

500, 30, and 0 km 

temperature profile 

U.S. Standard 

Atmosphere 1976 

sea-level pressure 

1 atm 

wavelength 

5 ym 
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All the data sets analyzed by the method of section 4 
contained 125 absorber amount values obtained using Eq. (64) 
and corresponding to 125 different Z Q values. For A 0 = 500 km 
the apparent zenith angles were chosen to have a constant 
angular spacing of 0.00879°, corresponding to a tangent height 
spacing less than 500 m over the entire range (0-50 km) . The 
distribution of the apparent zenith angles given in Table 2 
for the 0 km observer data sets were chosen to bias the data 
sets in favor of large apparent zenith angles to enhance the 
information content. This choice may affect the proportion- 
ality constants obtained for Eqs. (36) and (41). In general, 
the results obtained in this chapter are expected to be better 
than those obtained from the analysis of real data with an 
equivalent SNR. 

The Z 0 distribution for the 30 km observer contained 
60 values greater than 90° and these values were arranged so 
that the tangent height spacing was less than 500 m. In order 
to achieve this spacing the Z Q values were adjusted "by hand." 
Thus the apparent zenith angle spacing varied over the entire 
range above 90°. Below 85° an apparent zenith angle spacing 
of 5° was used. In the range 85°<Z o _<90 o the apparent zenith 
angle spacing was 0.10417°. 

The values assumed for the atmospheric parameters 
needed to construct the model atmosphere of section 4 are 
listed in Table 3. These were all assumed to be known 
exactly. 
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The effects of noise were simulated by adding normally- 

distributed random numbers with a mean of zero and a standard 

deviation o Q to the calculated absorber amounts. Since 

\ 

negative absorber amounts are unphysical, this sum was set 
equal to zero when the sum of the absorber amount and the 
random noise was negative. This preserved the total number 
N 0 of a m (Z 0 i) values upon which the vertical resolution of 
the retrieval depends. The SNR of each data set thus con- 
structed was defined as a /a where a was the largest 

max o max ^ 

absorber amount value in each data set. These synthetic 
data sets were used as "experimental" data and mixing-ratio 
values were retrieved by inverting Eq . (70) and constraining 

the mixing-ratio value above 50 km to zero. 

2.5.2 500 km observer 

Typical absorber amount data sets with an SNR of 100 
obtained from the constant mixing-ratio profile for the 
500 km and 0 km observers are shown in Figures 13 and 14 
respectively. A retrieval done with the data in Figure 13 
is shown in Figure 15. 

The percent uncertainty for the kth retrieved mixing- 
ratio value can be defined as 

100a k 

U k = , (73) 

K m ' 
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Total Absorber Amount (Arbitrary Units) 



□ - 


0 
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Zenith Angle (Decimal Degrees) 


Figure 13: Total absorber amount data with SNR of 100, 

calculated for the constant mixing-ratio profile 
for an observer at 500 km, as a function of the 
apparent zenith angle of arrival of the rays. 
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Total Absorber Amount (Arbitrary Units) 



Zenith Angle (Decimal Degrees) 


Figure 14: Total absorber amount data with SNR of 100, 

calculated for the constant mixing-ratio profile 
for an observer at 0 km, as a function of the 
apparent zenith angle of arrival of the rays. 
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Mixing Ratio 


Figure 15: A retrieval of the constant mixing-ratio profile 

performed on the data shown in Figure 13 for 
50 equally-spaced local mixing-ratio values. 
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where m is the value of the assumed constant mixing-ratio 
profile used to calculate the synthetic data. It was found 
that these uncertainty values could be decreased by reducing 
the number of values N R to be retrieved or by increasing 
the SNR or the number of data values N 0 . 

Retrievals of the constant mixing-ratio profile were 
performed for an observer at 500 km and different values of 
N R/ SNR, and N Q . As shown, in Table 4 the percent uncertainty, 
as defined by Eq. (73) , is approximately inversely propor- 
tional to the SNR. The data sets used to obtain these results 
were constructed from identical sets of random numbers scaled 
to give the SNR listed. 

As shown in Table 5, the percent uncertainty is inversely 

p 

proportional to N J . The results obtained for a data set of 
500 rays in Table 5 were obtained by combining four data sets 
each containing 125 values with an SNR of 100 but different 
added random noise. For each tangent height there were four 
absorber amount data values. Thus the information content 
of the data set was increased without affecting the vertical 
resolution of the retrieval. 

The dependence of the uncertainty on the number of 
equally-spaced mixing-ratio values Nr to be retrieved was 
determined by varying N R . The uncertainty values retrieved 
for N r values of 5, 10, 25, 50, and 100 are' shown in Figure 
16 as a function of altitude. The physically unrealistically 
large uncertainty values shown in Figure 16 were given to 
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Table 4: Some ratios of retrieved mixing-ratio uncertainties 

for N r = 100, N 0 = 125, A 0 = 500 km, and various 
SNR values . 




- 





SNR 


A (km) 

500 

100 

50 

10 

500 

500 

500 

500 

0 

1.0 

3.8 

7.2 

34 

10 

1.0 

4.0 

7.4 

36 

20 

1.0 

3.9 

7.5 

36 

30 

1.0 

4.0 

7.6 

36 

40 

1.0 

4.0 

7.6 

36 
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Table 5 : Some retrieved mixing-ratio uncertainties for 

Nr = 100, A 0 = 500 km, SNR = 100, and N 0 values 
of 125 and 500. 


A ( km ) 

N 0 = 125 

N 0 = 500 

U(N 0 = 125) 
U(N 0 = 500) 

0 

3.8% 

1.8% 

2.1 

10 

19 

9.2 

2.1 

20 

63 

31 

2.0 

30 

400 

190 

2.1 

40 

1800 

880 

2.0 
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demonstrate the quality of the fit to the straight lines 
over the entire altitude range. The straight lines in 
Figure 16 show that the uncertainty at any altitude is 
proportional to N^, since doubling N R doubles the uncertainty 
(note the logarithmic scale of the uncertainty axis) . 

In summary these results show that the uncertainty at 
any altitude obeys 


SNR N q 2 

The proportionality constant can be determined by examination 
of Figure 16. Since the uncertainty axis in this plot is 
logarithmic, the uncertainty increases exponentially with 
altitude. Noting that the slope of the lines in Figure 16 
is approximately 0.15 km - -*-, Eq. (74) becomes 


U (A, 500) 


50 


N R exp (0 . 15A) 


SNR N r 


(75) 


where 500 designates the observer altitude in kilometers. 

The scale height for the U.S. Standard Atmosphere 1976 
(28) is shown in Figure 17 as a function of altitude. The 
reciprocal of the mean value of H g is 0.14 km - -* - and the 
factor 0.15 km - -* - in Eq. (75) is thus nearly H s “l. For an 
isothermal atmosphere with constant molecular weight the 
atmospheric pressure at altitude A can be approximated by 
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U (Percent Uncertainty) 


Figure 16 : Uncertainty in the retrieved mixing-ratio values 

for an observer at 500 km obtained from data with 
SNR = 100, N 0 = 125, and N R = 100, 50, 25, 10, 
and 5. The straight lines are the corresponding 
functions from Eq. (2.75). 



SCALE HEIGHT (KM) 


Figure 17: The scale height of the U.S. Standard Atmosphere 

1976 (27) as a function of altitude. The vertical 
line represents the average scale height value of 
7.08 km . 





P (A) = P (o) exp (-A/H s ) 


(76) 


where P(o) is the pressure at sea-level. Thus the uncertainty 
in the mixing-ratio value at height A in terms of P(A) is 
approximated by 


[ 50 N r P (o) ] 

U (A , 500 ) = r 

.[SNR N q 2 P (A) ] 


(77) 


and it is inversely proportional to the pressure. Therefore 
the uncertainties in retrieved mixing-ratio profiles are 
primarily due to the rapid decrease in the air density of 
the hydrostatic atmosphere with altitude. If Eq . (77) is 
written in terms of the scale height 


50 N R P (o) exp ( A/H s ) 

U (A, 500) = £ . (78) 

SNR N q 2 

This formulation is in agreement with the result obtained 
from the Backus-Gilbert theory (Eq. (36)). 


2.5.3 0_ km observer 

For apparent zenith angles less than 90° the atmosphere 
is nearly plane-parallel. In this case 


a ij tanZ oj " a i , j+l tan z o j+1 


(79) 
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Thus for the observer within the atmosphere for absorber 
amounts corresponding to rays with apparent zenith angles 
less than 90° the aj_j elements for the jth ray will be near 
multiples of the elements of the coefficient vectors of all 
other rays with zenith angle less than 90°. For these 
conditions the {aj_j} matrix is nearly singular and the inver- 
sion of the matrix becomes extremely difficult. This near 
singularity could be reduced by reducing the number N R of 
parameters until the {aj_j} matrix is invertible. The inver- 
sion program (P9R) , beginning with the least important para- 
meter (determined by a calculated tolerance factor for each 
parameter) , eliminated parameters (starting with the smallest 
tolerance) until the (a^j) matrix is invertible. The retrieved 
mixing-ratio values obtained by this process are highly 
correlated due to the scaling in Eq. (79) . A more detailed 
explanation of the tolerance determination can be found in the 
BMDP manual (12). We have used a tolerance of 0.0001 in this 
work. For the 0 and 30 km observers the resultant number of 
retrievalbe parameters was less than N R except for N R less 
than five. 

Similar data sets to that shown in Figure 14 were analyzed 
for an instrument at 0 km. The program P9R found that it was 
not possible to obtain more than 10 local, mixing-ratio 
values in the 0 to 50 km region for all values of the SNR 
used (including the case of no added noise) . Although N max 
is essentially independent of the set of normally distributed 
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random numbers used to create the synthetic absorber amount 
data, the altitudes of the retrieved local, mixing-ratio 
values varied for each. data set analyzed. This was inter- 
preted as being due to the poor vertical resolution and the 
high correlation of the retrieved values. An example is 
shown in Table 6 for N 0 of 125, Nr of 50, and an SNR of 100 
for a constant mixing-ratio value of 1.0. The mixing-ratio 
values were retrieved from identical absorber amount data 
with different sets of noise. The vertical resolution was 
so poor above 3 km that the altitude variation of the 
retrieved values was more than 15 km for different experimen- 
tal noise but identical SNR. 

The retrieved standard deviations were also found to 
depend on the random numbers used to model experimental 
noise. The mean value of the standard deviation at some 
altitude A was found by averaging the retrieved standard 
deviation values from retrievals using 16 data sets con- 
structed in the same way as those for Table 6. As seen in 
Table 6, above 3 km the altitudes at which the mixing-ratio 
values were retrieved depended on the data set analyzed. 

Thus, for example, only once in sixteen retrievals was 
a mixing-ratio value retrieved at 40 km, while at 0 km the 
mixing-ratio value was always retrieved. At 0 km the average 
standard deviation could be calculated from 16 values but 
at 40 km the average consisted of just one value. These 
numbers of values used to obtain the averages are listed in 
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Table 6: Sample retrievals of constant mixing-ratio profile 

(m = 1.0) for Nq = 125, SNR = 100, Np = 50, and 
an observer altitude of 0 km. 
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Table 7. The existence of six retrieved values at 49 km 
was due to the constraint of the mixing-ratio to zero at 
.50 km. These values contain little information about the 
true profile. The spread of the standard deviations at any 
altitude about the average increases rapidly with altitude. 
Since the number of values retrieved decreases rapidly with 
altitude we have restricted our analyses to the region below 
30 km (below 15 km for an N R of 50) . The retrieval method 

used obtained a unique set of mixing-ratio and standard 

r 

deviation values determined mathematically by linear inver- 
sion with no regard to the possible physical meaning of the 
retrieved values. Thus, as shown in Table 6, physically- 
meaningless, negative mixing-ratio values were sometimes 
retrieved . 

The average uncertainty was obtained by using Eq . (73) 
where m is the value of the constant mixing-ratio used to 
calculate the absorber amount data. The average uncertainties 
for an N 0 of 125, an SNR of 100, and N R of 5, 10, and 50 are 
shown in Figure 18. The straight line estimates of the 
uncertainty values were obtained from 

U ( A, 0 ) = 45 U (A, 500) . (80) 

Shown in Table 8 are some of the average uncertainty values 
obtained for different values of the SNR. Comparison of the 
retrieved values in Table 8 shows that the uncertainty is 
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Table 7 : 


The altitude distribution of the mixing-ratio 
values retrieved from the 16 absorber amount 
data sets corresponding to an observer altitude 
of 0 km . 


A (km) 

Number of times retrieved in 
16 trials 

0 

16 

1 

14 

3 

12 

9 

3 

12 

2 

13 

2 

14 

1 

16 

o 

17 

2 

21 

1 

24-27 

4 

30 

1 

49 

6 
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U (Percent Uncertainty) 


Figure 18 : Average uncertainty in the retrieved mixing-ratio 

values for an observer at 0 km from retrievals 
from 16 data sets with SNR = 100, N p = 125, and 
N R = 50, 10, and 5. The straight lines are the 
corresponding functions from Eq. (2.80). 



Table 8: Some retrieved mixing-ratio uncertainties for 

N R = 50, N 0 = 125, A q =. 0 km, and various SNR 
values . 


A (km) 

SNR 


100 

1000 

5000 • 

0 

47% 

4.1% 

1.2% 

1 

97 

00 

• 

o 

2.4 

3 

250 

16 

4.9 

9 

590 



18 



100 

19 


400 


26 

4200 



23 



400 

29 


1700 
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approximately inversely proportional to the SNR. In Table 9 
are listed the values obtained for an N Q of 125 and 500, 
which indicate that the uncertainty is nearly inversely 
proportional to the square root of N Q . Thus the dependence 
of Eq. (80) on N Q and the SNR is essentially identical to 
that obtained for the 500 km observer. 

2.5.4 30^ km observer 

For an instrument at 30 km attempts were made to 
retrieve a total number Nr equally-spaced mixing-ratio values. 
It was found that the maximum number of mixing-ratio values 
retrievable below the observer was approximately equal to 
the number of rays with Z o >90° in the data set, but that 
very few values were obtained above the observer. 

The average uncertainties in the retrieved mixing-ratio 
values were calculated from 16 retrievals as described 
previously. Equation (80) also describes these retrieved 
uncertainties reasonably well; but at these altitudes there 
is a large scatter in the retrieved values. 

The uncertainties in the retrieved mixing-ratio profile 
below the 30 km observer show the same dependence of Nr, N q , 
and the SNR as for a 500 km observer. The values of U(A,30) 
for the atmosphere below 30 km are shown in Figure 19 for N Q 
of 125, an SNR of 100, and Nr values of 5, 10, 25, and 50. 

The straight lines calculated from Eq. (75), in this figure, 
show that the same expression can be used to describe the 
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Table 9 : Some retrieved mixing-ratio uncertainties for 

Nr = 50, SNR = 100, A Q = 0 km, and N Q values of 
125 and 500. 


A (km) 

Nq = 125 

N 0 = 500 

0 

47% 

20% 

1 

97 

39 

3 

250 

71 

9 

590 


19 


950 

26 

4200 


49 


65000 




Figure 19: Average uncertainty in the retrieved mixing-ratio 

values for an observer at 30 k.m obtained from 
retrievals from 16 data sets with SNR = 100, 

Nq = 125, and N R =50, 25, 10, and 5. The straight 
lines are the corresponding functions from 
Eg. (2.75). 
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uncertainties in the mixing-ratio profiles obtained at 
both 30 and 500 km. To illustrate this point some retrieved 
uncertainty values below the observer for various SNR 
values are listed in Table 10. In Table 11 a comparison 
is given of the uncertainties obtained for two different N 0 
values (all other parameters the same) . These show the 
same dependence as given in Eq. (80) . Therefore it is con- 
cluded that, Eqs . (75) and (80) describe the uncertainties 

in the mixing-ratio values below and above the instrument for 
all observing heights, for the case of an absorber with 
a constant mixing-ratio. 

It was found that, for other mixing-ratio profiles, the 
variations of a m j_(Z 0 ^) can be substantially different from 
those shown in Figures 13 and 14. This altered a . the 
effective SNR, and the values of the constants in Eqs. (75) 
and (80). However, the functional, form was qualitatively 
similar (i.e,, the uncertainty increased with height and 
depended on the SNR, Nr, and N Q in the same fashion as for 
the constant profile) . These other profiles are discussed 
later in this chapter. 

2.5.5- Summary 

In this long section we have determined the uncertainty 
in the retrieved mixing-ratio values as a function of alti- 
tude for the case of absorber amount data obtained for an 
absorber with a vertically constant mixing-ratio profile. 
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Table 10: Some retrieved mixing-ratio uncertainties below 

the observer for Nr = 100, N 0 .= 125, Aq = 30 km, 
and various SNR values. 


A (km) 

- 

SNR 


500 

100 

50 

10 

0 

0.80% 

3.9% 

7.3% 

35% 

5 

1.6 

6.3 

12 

56 

10 

2.6 

11 

20 

94 

15 

8.7 

35 

66 

310 

20 

13 . 

50 

95 

450 

25 

28 

110 

210 

1000 
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We have found that the uncertainty in these retrieved values 
increases exponentially with altitude, depends on the SNR of 
the data, the number of data values analyzed, and the number 
of mixing-ratio values retrieved, it is much larger above the 
observer than below, and it is independent of the observer 
altitude. Above the observer small uncertainties are only 
possible for small Nr, or large SNR, or large N 0 . For the 
region below the observer the uncertainty in the retrieved 
mixing-ratio values above 30 km was larger than 10% no 
matter what retrieval value spacing was used. Thus reasonable 
uncertainties in the retrieved mixing-ratio values above 
30 km also require small N R , or larger SNR, or large N 0 . 

These results were obtained for the case of a mixing- 
ratio profile retrieval from absorber amount data (equivalent 
width data) . These data represent only one piece of informa- 
tion from each spectrum. More information is contained in 
the shapes of the spectral lines. The uncertainty in re- 
trieved mixing-ratio values obtained from spectral curve- 
fitting is considered in the fourth chapter. 

The uncertainties in the retrieved mixing-ratio 
values could be reduced if a number of different spectral 
lines are analyzed.. If, for example, one spectral line with 
intensity S is used to determine the absorber amount for 
occultations corresponding to tangent heights less than 25 km 
and a second line with intensity lOOxS is used to determine 
the absorber amount for occultations corresponding to tangent 
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heights between 25 and 50 km the effective noise contained in 
the absorber amount data for tangent heights above 25 km will 
be less than if only the one line with intensity S were used 
for the entire region. Thus the overall effective SNR of the 
data will be reduced so that the uncertainty will also be 
smaller . 


2.6 THE VERTICAL RESOLUTION 
2.6.1 Resolution determination 

If N r equally-spaced, mixing-ratio values are retrieved 
over an altitude interval AA (=50 km) the corresponding 
vertical resolution can be defined as 

Res = ^ . ' (81) 

n R 

It has been found that if the SNR is sufficiently high. Nr 
can be increased to nearly N 0 and the ultimate resolution in 
a mixing-ratio profile retrieval is 


r - 2AA 
ULT . N c 


(82) 


Equations (75) and (80) and the straight lines in Figures 
16, 18, and 19 show the dependence of the uncertainty of the 
retrieved mixing-ratio values on altitude for a particular 
vertical resolution defined by the spacing of the retrieved 
mixing-ratio values. But in order to determine the optimum 
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distribution of retrieval values in the analysis of a given 
set of a data, it is more useful to know the dependence of 
the vertical resolution on altitude for a fixed retrieval 
value uncertainty at all altitudes. 

By combining Eqs . (75) and- (81) the profile resolution 

Res B (A) below the observing instrument is 

Res B (A) = IQOA A exp (0 . 15A) # (83) 

No U SNR 

Some examples are shown in Figure 20 for R ULT = 1.0 km and 
constant values of U of 0.1, 1, 5, and 25%. A vertical reso- 
lution greater than 25 km is unrealistic for the altitude 
range of 50 km sampled. Figure 20 shows a vertical resolution 
up to 100 km in order to describe the dependence of Res B on 
altitude. The cut-off at 1 km is due to the limit r ult = 1 km. 
The form of Eq. (83) is in agreement with Eq. (41) obtained 
from the Backus-Gilbert theory. The vertical resolution is 
inversely proportional to the retrieval value uncertainty. 

The corresponding relation for the atmosphere above the 
observer, obtained from Eq. (80) is 

Res A (A) = 4 5Res B (A) . (84) 

Examples of the dependence of the resolution Res A (A) on 
height for the case of R^ =1.0 km and an uncertainty CJ of 
5, 25, and 100% are shown in Figure 21. Even for the 
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R (Vertical Resolution KM) 


Figure 20: Vertical resolution of the retrieved profiles 

obtained from Eq . (2.83) for the region below 
the observer for SNR = 100, R ULT = 1.0 km, and 
U = 25.0, 5.0, 1.0, and 0.1%. 




U = IOO%> 


11=25*5 


U=5%; 


I 10 100 

R (Vertical Resolution KM) 


Figure 21 : 


Vertical resolution of the retrieved profiles 
obtained from Eq. (2.84) for the region above 
the observer with SNR = 100, R ULT = 1.0 km, 
and U = 100, 25, and 5%. 


optimistically large values of 125 and 100 for N 0 and the 
SNR respectively, it is not possible to obtain significant 
resolution above the observer if an uncertainty of 5% or 
better in the retrieved mixing-ratio values is desired. 

Indeed, to obtain a vertical resolution of 10 km at a height 
of 20 km above the instrument, together with an uncertainty 
of 5% in the retrieved mixing-ratio values, an SNR of 10 4 is 
required for the case of a constant mixing-ratio retrieval. 

In summary, if no a priori information about the vertical 
mixing-ratio profile is assumed, the quality of the retrieved 
profiles degrades nearly exponentially with altitude, is 
considerably better for the atmosphere below than above the 
observer, and is independent of the instrument height. 

2.6.2 The spacing of retrieval parameters 

The previous results suggest that, by spacing the local, 
mixing-ratio values to be retrieved according to the expected 
vertical resolution and desired uncertainty, it is possible 
to reduce significantly the number of atmospheric layers used 
in the calculations and the number of local, mixing-ratio 
values to be retrieved, and hence the amount of computer time 
required for the analysis. 

The vertical resolution defined in Eq. (81) was given as 
twice the mixing-ratio value spacing in the retrieval. 
Similarly, the spacing of adjacent parameters at altitude A 
is given by one half the vertical resolution at altitude A 
f rom Eq . (83). 
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This suggests that if the first mixing-ratio value is 
to be retrieved at = 0 km, the second mixing-ratio value 
should be retrieved at approximately A 2 = R 0 /2 where, from 
Eq. (83) 


Ro = 


100AA 
U SNR 


(85) 


in order to give a constant uncertainty U in the retrieval 
values. Similarly, the altitude of the third retrieval 
value is 


A 3 = A 2 + exp (0 .15A2) 


( 86 ) 


and the altitude of the kth mixing-ratio retrieval value 
should be approximately 


A k“ A k -1 



exp ( 0 . 15A k-1 ) 


(87) 


In principle, Eq. (87) gives the altitude of the kth 
retrieval parameter in terms of the vertical resolution at 
R q . Equation (87) cannot be expressed in a simple closed 
form, nor it is well-behaved when approximated by the Taylor 
series . 

It can be used iteratively to determine the spacing of 
retrieval values in order to give a constant uncertainty U. 
The appropriate retrieval parameter altitudes for A 0 = 500 km. 
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AA=50 km, N o = 100, U=l%, and an SNR of 10 0 and 500 are 

shown in Table 12. A retrieval performed on data obtained 

from a constant mixing-ratio profile, for the conditions of 

Table 12 and an SNR of 100 are given in Table 13 . The 

uncertainties in the retrieved values at higher altitudes 

are substantially reduced for the variable spacing retrieval 

and are approximately independent of altitude. A comparison 

of the uncertainties in this retrieval and those obtained for 

a constant spacing show that the appropriate distribution of 

retrieval values allows the retrieval of mixing-ratio values 

/ 

over the entire region of interest with nearly constant 
uncertainty. 

The distribution of mixing-ratio values for a 0 km 
observer is substantially worse. For A A = 50 km, N o = 100, 
and U = 1% an SNR of 225 is required to retrieve two mixing- 
ratio values, and an SNR of 1200 is necessary to retrieve 
three values. Retrievals performed by Zander et al (29) on 
absorption spectra with an SNR of 300 obtained at Kitt Peak 
in order to determine the column abundances of FC-11, FC-12, 
and FC-22 show that not more than one well-determined mixing- 
ratio value is obtainable. Thus it is concluded that our 
results are in general agreement with experiment. 

A parameter spacing which varies with altitude has been 
used by Goldman et al (30) to retrieve mixing-ratio values 
from spectra taken at ground-level. In this case the uncer- 
tainties were reduced by using a scheme similar to the one 
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Table 12: Sample retrieval profiles which give an uncertainty 

in the mixing-ratio values of 1% at all altitudes 
for AA = 50 km, N 0 = 100, and an SNR of 100 and 
500. 
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Table 13: A sample mixing-ratio profile retrieval of the 

constant profile of m = 1.0 using the retrieval 
profile of Table 12 for SNR = 100, N 0 = 100, 

A A = 50 km. 


Retrieved uncertainty 
values 

Uncertainty values 
predicted from Eq. (54) 
for Nr = 5 

A (km) 

m k 

Uk 

A (km) 


Uk 

0 

1.03 

2.6 

0 


0.2 

2.5 

0.99 

1.7 

10 


1.1 

6.0 

0.99 

1.5 

20 


5.0 

12.5 

1.02 

2 .0 

30 


22 

28.5 

0.98 

6.0 

40 


101 
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described here. The results obtained in this work suggest 
that the representation of mixing-ratio retrievals can be 
improved by suitable spacing of the mixing-ratio values. 

2.7 RETRIEVALS OF OTHER PROFILES 

The nature of the trade-off, due to the experimental 
noise in absorber amount data, between the uncertainty in 
the retrieved mixing-ratio values and the vertical resolution 
of other assumed mixing-ratio profiles has been examined by 
considering a number of profiles. The standard values of 
Table 3 were used with an observer altitude of 500 km. As 
was the case for the previous retrievals no a priori assump- 
tions were made about the shape of the profile to be retrieved. 

Two of the profiles considered, the quadratic profile 
of Figure 22 and the quadratic spike shown in Figure 23, were 
chosen for study since many atmospheric molecules have similar 
mixing-ratio profiles (31) . The quadratic profile is similar 
to the profile of HNO 3 and the quadratic spike similar to the 
profile of NO 2 . 

In addition to these smoothly varying profiles, we have 
constructed a discontinous profile in order to test the 
uncertainty-resolution trade-off found previously. This 
profile, shown in Figure 24, consisted of six 1 km-thick 
segments with constant mixing-ratio values. In between 
these segments the mixing-ratio values were zero. Thus all 
the absorber was concentrated in six discrete regions of 
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50 



Figure 22: A mixing-ratio profile which varies quadratically 

with altitude and is a maximum at 25 km. 
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HEIGHT (KM) 



Figure 23 : A mixing-ratio profile which varies quadratically 

with altitude between 15 and 40 km and is a 
maximum at 37.5 km. 
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the atmosphere separated by 3 km in altitude (a highly 
unphysical situation) . The mixing-ratio values of each 
segment were chosen such that a vertical column in each 
segment contained the same number of absorber molecules. 

The data sets constructed from these mixing-ratio 
profiles contained absorber amount values corresponding to 
125 apparent zenith angles of constant spacing with tangent 
heights between 0 and 50 km (the range of tangent heights 
for the "comb" profile of Figure 24 was 15 to 40 km) . The 
tangent height spacing was less than 500 m for all tangent 
heights (for the "comb" profile the spacing was less than 
250 m) . The observer altitude was 500 km. 

The data set calculated using the quadratic profile is 
shown in Figure 25. Random noise was added to the calculated 
absorber amounts to give an SNR of 200. The data set obtained 
from the quadratic spike with an SNR of 100 is shown in 
Figure 26. 

The absorber amount data shown in Figure 27 calculated 
from the comb profile indicate the existence of the discrete 
non-zero layers. The peaks allow the separation of the 
layers if the noise is sufficiently small. The SNR of the 
data in Figure 27 is 50. 

One hundred local, mixing-ratio values with a spacing 
of 500 m between 0 and 50 km were retrieved from the data 
shown in Figure 25. The retrieved profile is shown in Figure 
28. As expected the uncertainty in the retrieved values 
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TOTAL ABSORBER AMOUNT (ARBITRARY UNITS) 



Figure 25: The total absorber amount as a function of the 

tangent height of the ray with an SNR of 200 
obtained from the quadratic mixing-ratio profile 
and an observer altitude of 500 km. 
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Figure 27: The total absorber amount as a function of the 

tangent height of the ray with an SNR of 50 
obtained from the comb mixing-ratio profile and 
an observer altitude of 500 km. 
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MIXING RATIO (ARBITRARY UNITS) 


Figure 28: A retrieval obtained from the data of Figure 25 

for a mixing-ratio spacing of 0.5 km. 
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increases rapidly with height. Above 35 km it is so large 
that no useful values are obtained. A retrieval performed 
on these data for a mixing-ratio value spacing of 5 km is 
shown in Figure 29. The uncertainty is greatly reduced and 
the altitude range of useful values is extended to 45 km. 

These profiles show that the uncertainty in the mixing-ratio 
values , can be reduced by increasing the spacing between para- 
meters in agreement with the results obtained for the con- 
stant profile in Eq. (75) . The effect of reducing the number 
of parameters is to increase the amount of data per parameter 
retrieved. The retrievals shown in Figures 28 and 29 show 
no systematic errors in the retrieval. At low altitudes the 
retrieval fits the true shape extremely well. At higher 
altitudes departure from the true profile are due only to 
the noise in the synthetic data. No obvious calculational 
errors were found in any of the retrievals done for this work. 

Similar results were found from analysis of the quadratic 
spike data shown in Figure 26. The uncertainty increased 
rapidly with altitude but could be reduced by increasing 
the parameter spacing. Retrievals obtained for a spacing of 
500 m are shown in Figure 30 and those obtained for a 5 km 
spacing are shown in Figure 31. 

Ten data sets similar to those shown in Figure 26, 
obtained from the comb profile, were constructed by adding 
different sets of random noise so that in all cases the SNR 
was 50. 
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Figure 29: A retrieval obtained from the data of Figure 25 

for a mixing-ratio spacing of 5 km. 
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Figure 30: A retrieval obtained from the data of Figure 26 

for a mixing-ratio spacing of 0.5 km. 
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Figure 31: A retrieval obtained from the data of Figure 26 

for a mixing-ratio spacing of 5 km. 
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Retrievals were made by using these ten data sets and 
a parameter spacing of 250 m. At each altitude the average 
parameter value together with its standard deviation was 
calculated from the ten retrievals. The effect of averaging 
the ten sets of mixing-ratio values with SNR of 50 was 
identical to one retrieval made on data with SNR of about 
160 as expected. The average mixing-ratio values are shown 
in Figure 32 together with the standard deviations. The 
retrieved values are seen to be accurate at all altitudes. 

The standard deviations, shown in Figure 33, increase with 
altitude (except above 40 km, where the standard deviations 
are smaller than expected due to the constraint on the mixing- 
ratio value at 50 km) as was the case for the constant 
mixing-ratio profile retrievals. 

Similar averages were obtained using the same data sets 
for a parameter spacing of 2.5 km. The calculated averages 
and standard deviations are shown in Figures 34 and 35. 
Comparison of Figure 33 and 35 show that although the uncer- 
tainty in the retrieved values can be reduced by increasing 
the parameter spacing, it is accompanied by a degradation in 
the vertical resolution. 

The previous analysis suggests that if no a priori 
assumptions are made about the shape of the profile, the 
uncertainty in the retrieved mixing-ratio values is given by 
Eq. (75) or Eq. (80) and the vertical resolution of the 
retrieved profile is given by Eq. (83) or Eq. (84) regardless 
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Figure 32: The average of ten retrievals obtained from ten 

data sets similar to that shown in Figure 27 for 
a mixing-ratio spacing of 0.5 km. 
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of the shape of the profile to be retrieved. For a particular 
profile the constants in these equations will have different 
values than those defined for Eqs . (75), (80), (83), and (84) 
but the functional form of these equations remains the same. 
These constants can be determined for any profile of interest 
by the methods already described. 

Analysis of data sets obtained for various profiles and 
an observer altitude of 0 km show, as expected from Eq. (84) , 
that there is very little ability to discriminate a mixing- 
ratio profile above the observer. An example shown in Figure 
36 was obtained from data calculated using the quadratic 
profile (with no added noise) with an N R of 50 and an N 0 of 
125. As seen in this figure no useful information is 
retrieved above 3 km and the vertical resolution is very 
poor. The number N R of parameters retrieved was found to be 
independent of the profile studied and in all cases the 
uncertainties were large and increased rapidly with altitude 
and the vertical resolution was very poor. 


2.8 SUMMARY 

In this chapter we have obtained expressions for the 
uncertainty in retrieved mixing-ratio profiles, if no 
assumptions are made about the shape of the profile to be 
retrieved, both above and below the observer altitude. These 
uncertainty expressions suggest that many investigators may 
have underestimated the errors in their reported retrievals 
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from solar spectra. By the methods outlined in this chapter 
it is now possible to estimate the uncertainty in a mixing- 
ratio profile retrieval based entirely, on the experimental 
conditions and not on any a priori assumptions about the 
errors. 

The uncertainty in retrieved mixing-ratio values was 
found to increase exponentially with altitude, to decrease 
with increasing SNR, to increase if more values were re- 
trieved, to decrease if more data values were analyzed, it 
was much larger above the observer than below, and it was 
independent of the observer altitude. These relations were 
found to be true for all the mixing-ratio profiles studied. 
Thus the uncertainty expressions (with slight modification) 
can be used to express the retrieval of any mixing-ratio 
profile from equivalent width data. 

The vertical resolution of retrieved mixing-ratio 
profiles degraded exponentially with altitude and was much 
worse above than below the observer. The vertical resolution 
expressions can be used to determine the appropriate spacing 
of the retrieval values in order to obtain an uncertainty in 
the retrieved values which Is constant with altitude. 

Thus the effects of experimental noise on the analysis 
of equivalent width data obtained from a single spectral 
line can be described by the equations obtained in this 
chapter. Similar analyses can also be used to describe 
other experimental situations , such as the analysis of 
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equivalent width data obtained from several spectral lines in 
each spectrum or equivalent width data obtained from spectral 
lines with different intensities (each intensity is used to 
obtain equivalent widths with corresponding tangent heights 
in a particular region of the atmosphere) . The additional 

information contained in the spectral line shape is considered 

* 

in the fourth chapter. 

In the next chapter we investigate the effects of syste- 
matic errors in the describing- physical parameters on the 
retrieval of mixing-ratio profiles. 
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Chapter III 
SYSTEMATIC ERRORS 


3.1 INTRODUCTION 

In the previous chapter we determined functional 
relationships for the uncertainty and the vertical resolution 
of mixing-ratio profiles retrieved from analysis of synthe- 
tic absorber amount data. In order to do this we assumed 
that the thermal structure of the atmosphere and the orien- 
tation of the observing instrument were exactly known. 
Systematic errors in the parameters necessary to describe 
the atmospheric structure and the position of the observer 
will cause additional uncertainties in the retrieved mixing- 
ratio values. In this chapter we examine the effects of 
these systematic errors on the retrieval of mixing-ratio 
profiles from absorber amount data. 

Other investigators have examined the effects of errors 
in the surface pressure, temperature profile and spectral 
line intensities, the effects of errors due to neglect of 

I 

refraction and the Doppler shift, and the smearing of spectral 
data due to the finite field of view of the sun on the inter- 
pretation of satellite radiance measurements in order to 
determine mixing-ratio profiles of atmospheric constituents 
(1-3) and the atmospheric temperature and pressure profiles 
(4,5). The effect of error in the atmospheric temperature 
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profile has also been examined for the retrieval of the CO 
mixing-ratio profile from solar absorption spectra obtained 
from balloons (6) . No studies have been performed on the 
effects of systematic errors on the retrieval of mixing-ratio 
profiles (with no a priori assumptions about the profile 
shape) from analysis of equivalent width data. 

The types of systematic errors considered in this chapter 

were 

1) those due to incorrect physical modeling of the 
slant path problem, 

2) those due to incorrect assumed values for para- 
meters describing the thermal structure of the 
atmosphere, and 

3) those due to incorrect assumed vaiues for para- 
meters describing the orientation of the observing 
instrument with respect to the earth's surface and 
the source of the radiation. 

The systematic errors in the retrieved mixing-ratios consid- 
ered in each category were determined by the slant path 
calculating program (7,8) used for the optical airmass com- 
putation. The slant path program, as described in the previous 
chapter, required the specification of the temperature at 
discrete altitudes , > the pressure at some base level, the 
apparent zenith angle of the incoming radiation, the altitude 
of observation, the latitude of observation, and the wave- 
length of the observed radiation. Thus the Type 1 errors 
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consisted of spacings of the layer break points so large that 
the numerical integrations of Eqs . (A. 95) - (A. 97) were not 
accurate, and the effects of neglecting atmospheric refraction 
(included here were errors in the wavelength of the observed 
radiation since the index of refractions of air given by 
Edlen's formula (9) is wavelength dependent). The Type 2 
errors considered were in the assumed temperature profile and 
in the assumed sea-level pressure. Type 3 errors included the 
apparent zenith angle of observation, the altitude of obser- 
vation, and the latitude of observation. These errors are 
summarized in Table 14. 

As shown in Figure 37, a mixing-ratio profile can be 
obtained by comparing the absorber amount data along several 
ray paths to the corresponding optical airmasses along 
identical ray paths through the atmosphere. For the case of 
a constant mixing-ratio profile, as shown in Figure 37 b) , 
the mixing-ratio at any altitude can be found by taking the 
ratio of the absorber amount to the optical airmass for the 
ray with tangent height corresponding to the desired altitude.* 
If there is a systematic error in the defining parameters, 
the optical airmass values are incorrectly determined (see 
Figure 38) . Thus the ratio of the correct absorber amounts to 
the incorrect airmasses gives the incorrect retrieved 

*This is not the retrieval method used here. Rather each 
path is broken into N R segments and all the N 0 data are analyzed 
simultaneously. 
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Table 14: The systematic errors considered in this work. 


Type 1 

due to incorrect modelling of the physical problem 

layer thickness 
refraction 

Type 2 

due to incorrect values for defining atmospheric parameters 

P (0) 

T (A) 

Type 3 

due to incorrect values for defining observer 
geometry parameters 







b) 

5 



* 


a 

< 


m *¥ 


MIXING RATIO 


Figure 37 : Schematic retrieval of a constant mixing-ratio 

profile by comparison of the optical airmass to 
the absorber cunount. The mixing-ratio at altitude 
A in b) is obtained by the ratio of the two 
curves in a) at tangent height A. 
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AIRMASS 


ABSORBER AMOUNT 


TANGENT HEIGHT (KM) 


TRUE 

PROFILE 


RETRIEVED 

PROFILE 


MIXING RATIO 


Schematic retrieval of a constant mixing-ratio 
profile in the presence of a systematic error. 
The error causes a shift in the optical airmass 
curve a) which gives an erroneous profile in b) 


mixing-ratio profile of Figure 38 b) . It is instructive to 
estimate the error introduced into the retrieved mixing-ratio 
values by errors in the defining parameters and to estimate 
the critical altitudes where these errors become smaller than 
the uncertainty due to experimental noise. 

The effects of systematic errors were estimated by 
analyzing data sets consisting of 125 synthetic absorber 
amounts, corresponding to rays with 125 different apparent 
zenith angle values. These absorber amounts were calculated 
by using the standard values of Table 15 and a mixing-ratio 
profile which was constant with altitude. No noise was 
added to these data before the analysis. The aj_j values of 
Eq. .(2.70) were calculated assuming an error in each parameter 
of Table 15 independently. Mixing-ratio values were retrieved 
by analyzing the synthetic absorber amount data with the linear 
inversion program (10) described in the previous chapter. The 
method of analysis used in this work was a direct linear 
least-squares method which determined all the retrieved mixing- 
ratio values simultaneously. Unlike other methods, such as 
the "onion skin peeling" technique (11,12), this precludes the 
possibility that the effects of the systematic errors are due 
to errors introduced by the retrieval method. The results 
obtained agree qualitatively with those predicted by variation 
of Eq. (2.70) and the use of some simplifying assumptions. 

We have considered only the case of a 500 km observer, 
although the analyses suggest that similar errors can be 
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Table 15: Standard atmospheric values used for investigation 

of systematic errors . 


Parame ter 

Assumed value 

sea-level pressure 

1 atm 

temperature profile 

U.S. Standard 

Atmosphere 197 6-*-^ 

apparent zenith angles 

110° - 112° 

observer. altitude 

500 km 

latitude of observer 

45° 
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expected for observers at other altitudes. We have also 
briefly examined the effects of smearing of the retrieved 
information due to the motion of the observer (14) or the 
finite size of the solar disc viewed (15) . 

In the next section the effects of systematic errors on 
the calculation of the optical airmass encountered along the 
ray path are discussed. In the following section we consider 
the much more interesting problem of the effects of these 
errors on the retrieved mixing-ratio profiles. In the fourth 
section estimates are given for the combined error in a 
mixing-ratio profile and the critical altitude at which the 
systematic errors are comparable to the uncertainty due to 
noise. In the final section a method of internally calibrating 
the data in order to eliminate these errors is described. 

The relative merits of internal calibration as opposed to using 
externally determined parameters are discussed. 

3.2 THE EFFECTS OF SYSTEMATIC ERRORS ON THE AIRMASS 
3.2.1 Type 1_ - layer thickness 

The calculation time necessary for the analysis can be 
reduced by decreasing the number of atmospheric layers. This 
savings may result in a loss of accuracy in the calculation of 
the optical airmass traversed. 

The effects of the layer thickness on the optical airmass 
calculation were estimated by calculating optical airmasses 
for the standard values of Table 15 with layer thicknesses 
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ranging from 25 m to 20 km. Three apparent zenith angles were 
chosen, corresponding to rays with tangent heights of 48.5 km, 
12.0 km and 2.6 km. The percent error in the calculated total 
optical airmass traversed by the rays was defined as 


% error 


100 


a ~ a 25 \ 
i a 2 5 I 


( 1 ) 


where a 25 is the optical airmass calculated for a layer 
thickness of 25 m. 

These errors are shown in Figure 39 as a function of the 
layer thickness. For a tangent height of 2.6 km the error is 
insignificant for layer thicknesses up to 20 km. However, for 
a tangent height of 48.5 km a layer thickness larger than 600 m 
causes errors greater than 1%. 

The narrowness of the layers required to maintain accuracy 
for paths with large tangent heights was unexpected since the 
atmospheric density above 50 km is three orders of magnitude 
smaller than that near the earth's surface and the effects of 
refraction were expected to be negligible. The small layer 
thickness required for accuracy is believable due to the dis- 
continuities in the assumed temperature profile, shown in 
Figure 40, near the tangents of 12 and 48.5 km. The assumption 
that the . temperature gradient is constant in the tangent height 
layer is not true for these tangent heights for layer thick- 
nesses larger than 1 km. 
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ERROR 



LAYER THICKNESS (METERS) 


Figure 39: The error in the calculated optical airmass due 

to improper layer thickness for tangent heights 
of 2.64,- 12.04, and 48.48 km. 
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HEIGHT (KM) 



Figure 40: The U.S. Standard Atmosphere 1976. 
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The consistency of each of the curves in Figure 39 sug- 
gests that the round-off errors in the double-precision 
calculations were small. It also appears that, for a ray with 
a 2 km tangent height, only 5 layers between 0 and 100 km 
altitude are required to estimate the total optical airmass 
accurately although larger relative errors may occur in the 
airmass along specific segments of the path. For a ray with 
a 50 km tangent height, a layer spacing of 500 m near the 
tangent height is required if accuracies of 1% in the total 
airmass are needed. We have chosen a layer thickness of 250 m 
for the remainder of the calculations of this section. The 
problem of a layer thickness which varies with altitude has 
not been considered. 

For the analysis of absorber amount data containing 
experimental noise the results of the previous chapter suggest 
that, due to the inherent lack of information, contained in the 
data because of the noise, the layer thickness can often be 
increased from the values indicated by this analysis. However, 
the presence of the structural features obtained in the 
retrievals of a constant mixing-ratio profile (these retrievals 
are discussed in the next section) imply that even small layer 
thicknesses can introduce artifacts into the retrieved profile. 

3.2.2 Type 1 - refraction index of air 

The index of refraction of air is calculated from Edlen's 
formula and the Lorenz-Lorentz formula (16). For T = 288 K, 
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P = 1 atm and X = 5 vim, at sea-level 


(n G - 1) = 2.727 x 10 


( 2 ) 


Since (n-l)<<l at all altitudes, and from the Taylor expan- 
sion of the Lorenz-Lorentz formula, (n-l)ap, the importance 
of refraction decreases rapidly with altitude due to the 
exponential decrease in the air density. Thus it might be 
expected that refraction has little effect on the path tra- 
versed by solar radiation especially at high altitudes. 

In order to measure the effects of refraction, the air- 
mass was calculated for the standard values of Table 15 for 
rays with tangent heights between 0 and 90 km. These values 
were compared to optical airmasses calculated for the same 
values of the apparent zenith angle but (n o -l)=0. This 
causes the index of refraction at all altitudes to be identi- 
cally equal to one. Some of the total airmasses obtained for 
these conditions are compared in Table 16. The error was 
calculated for each tangent height by 


% error 


100 


a 

unre 

a 


f ~ a ref 
ref 


(3) 


where a unre ^ was the optical airmass assuming no refraction 
and a re f was the airmass if refraction was included. The error 
is shown in Figure 41 as a function of the tangent height of 
the refracted ray. It is less than 1% for rays with tangent 
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Table 16: A comparison of slant paths calculated including 

the effects of refraction to those neglecting 
refraction. 



Airmass 


Tangent Height (km) 

D 

with 

refraction 

with 

no refraction 


refraction 

no refraction 

110.0 

0.004 

0.0004 


85.75 

85.75 

110.2 

0.0014 

0.0014 


77.52 

77.52 

110.4 

0.0048 

0.0048 


69.20 

69.20 

110.6 

0.015 

0.015 


60.80 

60.80 

110.8 

0.043 

0.043 


52.33 

52. ,33 

111.0 

• 0.12 

0.12 


43.77 

43.78 

111.2 

0.42 

0.42 


35.13 

35.15 

111.4 

1.6 

1.6 


26.39 

26.44 

111.6 

6.5 

6.3 


17.46 

17.65 

111.8 

28 

24 


8.04 

8.79 
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% ERROR 


Figure 41: The error in the calculated optical airmass due- 

to neglect of the effects of refraction. 
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heights above 25 km but it increases rapidly as the tangent 
height decreases. The small difference in tangent height, 
between the refracted and unrefracted ray, causes a large 
change in the optical airmass. Refraction can be ignored 
above 25 km if accuracies of the order of 1% are sufficient. 
The results shown in Figure 41 agree qualitatively with those 
obtained by Snider and Goldman (8) for an observer at 40 km. 
Wavelength dependence 

Although the wavelength of spectral features can be 

C 

determined typically to 1 part in 10 the path traversed by 
radiation in the visible spectrum may be substantially dif- 
ferent than that traversed by infrared radiation, for rays 
arriving at the instrument at the same zenith angle, due to 
the wavelength dependence of the refractive index of air. 

If the observing instrument is aligned with the sun using the 
visible image of the sun, the apparent position of the infra- 
red image may be displaced from the visible image. 

The effects of changes in the wavelength were estimated 
by determining the variation of (n Q - 1) between 0.1 um and 
20 um as shown in Table 17. Over most of the region this 
variation is less than 1% with respect to the value at 5 um 
.although it increases rapidly at wavelengths below 0.5 um. 

The systematic errors in total airmass traversed due to 
changes in the wavelength given by the relation 


% error 


100 


a 


- a„ 
3Uin 

5um i 


(4) 
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Table 17: Some values of (n 0 - 1) of Edlen's formula for the 

reference T 0 = 288.16 K, P 0 = 1 atm, and various 
wavelengths . 




are shown in Figure 42 as a function of wavelength for three 
tangent heights (48.5 km, 12 km, and 2.6 km) . For rays with 
tangent heights less than 50 km, the errors were less than 1% 
for X > 0.4 pm. For the infrared region (X>0.7 pm) the error 
was less than 0.5% for all tangent heights less than 50 km. 

This implies that if the instrument is aligned with the sun in 
the visible, then there is little error in the estimated total 
optical airmass traversed by the infrared radiation by using 
the apparent zenith angle of visible wavelengths. 

3.2.3 Type 2_ - sea-level pressure 

The optical airmass encountered by a ray passing through 
the atmosphere is given by 

N 

a = l p 6 £ -1 , (5) 

j=l J J 

where Pj is the local air density in layer j, <55, j is the 
length of the ray path through layer j , and N is the number of 
layers. If the air density is written as some function f(A) of 
the altitude 

Pj = p(o)f(Aj) , (6) 

where f(Aj) is a function of temperature and altitude but not 
of the sea-level density, then Eq. (5) becomes 
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ERROR 


+ 2.64 

o 12.04 
□ 48.46 



Figure 42: The error in the calculated optical .airmass due 

to error in the radiation wavelength with respect 
to 5 ym for tangent heights of 2.64, 12.04, and 
48.46 km. 
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N 

a = Z p (o) f (A-; ) 5 1 • 

j=l 3 


(7) 


A difference in sea-level density Ap(0) results in a varia- 
tion of the optical airmass Aa given by 


N 

Aa = Z Ap (o)f (A.) 61 . , 


j=l 


3 J 


( 8 ) 


provided the temperature profile remains constant, or 


. , . N 

Aa = Z p (o) f (Aj) . 


(9) 


j=l 


Thus the fractional change in the optical airmass is 


Aa _ Ap (o) 
a • p(o) * 


( 10 ) 


If we assume that the pressure P(A) at altitude A is given by 


P (A) = P (o ) f (A) 


( 11 ) 


then 


Ap (o) _ AP (o) 

P (o) P (o) 


and 


Aa _ AP (o) 
a P(o) 
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Thus if the sea-level pressure is uncertain by x percent the 
total optical airmass is also in error by x percent. This, 
of course, does not show the resultant errors in a mixing- 
ratio profile retrieved using these erroneous airmasses. 

This problem is dealt with in the next section. 

The predictions of Eq. (13) were tested by comparing the 
optical airmass obtained from the standard values of. Table 15 
to those obtained for various errors in the base level pres- 
sure. The systematic error due to errors in the sea-level 
pressure was defined by 



where a s was the total airmass obtained from the standard 
values, and is shown in Figure 43 for four tangent heights. 

As expected, the error in the optical airmass is nearly 
independent of tangent height for a pressure error of less 
than 5%. The. small increase with decreasing tangent height 
from the predictions of Eq. (13) are due to the effects of 
refraction which were not included in Eq. (8) . There 
appears to be little advantage to moving the reference 
pressure base level up to the tangent height since the error 
in the calculated optical airmass is essentially independent 
of the tangent height. In fact, the pressure at the tangent 
height may be less certain than the corresponding ground- 
level pressure. 
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0 2.5 

APo (PERCENT OF P 0 ) 

Figure 43: The error in the calculated optical airmass due 

to error in the sea-level pressure for tangent 
heights of 2.64, 5.87, 12.04, and 48.48 km.. 
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For the long ray paths through the atmosphere encountered 
in occultation spectroscopy the sea-level pressure directly 
below these paths may vary appreciably along the ray path. 

As shown in Figure 44 the length of the grazing ray path 
through the lowest 50 km of the atmosphere is approximately 
1000 km. As seen in the surface weather map for the United 
States on April 9, 1980 (17) given in Figure 45, the typical 
variation of the sea-level pressure over such distances is 
typically of the order of 5%. Thus from Eq. (13) it is 
expected that the error in the airmass determination could be 
as large as 5%. 

Because of the rapidly decreasing air density with 
altitude most of the airmass is encountered in the vicinity 
of the tangent height. The error in the airmass determination 
due to base level pressure errors can be reduced by using 
the base level pressure directly below the tangent height. 

3.2.4 Type 2_ - temperature profile 
Horizontal and seasonal variations 

We have assumed that the thermal structure of the atmos- 
phere was correctly represented by the U.S. Standard Atmosphere 
1976 and that the structure was "horizontally " uniform along 
the ray path. As shown in Figure 44, and noted previously, 
the path length through the lowest 100 km of the atmosphere, 
for a grazing ray, is approximately 2000 km. This path 
length corresponds to a maximum latitude change of 10°. 
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< 2500 KM > 

«I000 KM> 



Figure 44: The lengths of grazing ray paths through the 

atmosphere . 
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Figure 45: Surface weather map for April 9, 1980 showing the 

sea-level pressure. 
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Even in the ideal case, where the atmosphere at some point 
along the path, preferably the tangent height, is adequately 
described by the U.S. Standard Atmosphere, departures from 
this temperature profile can occur along the paths of grazing 
rays . 

Table 18, obtained from Kantor (18), gives the mean 
monthly latitudinal density gradients for the mid-latitude 
atmosphere. These density gradients are due both to changes 
in sea-level pressure and the changing temperature profile. 

For a 2000 km path length the mean variation of the density 
at any altitude is about 5% although the instantaneous varia- 
tion can be much greater. 

Errors in the estimated airmass may also arise if an 
improper atmospheric temperature profile is used. The U.S. 
Standard Atmosphere, used in this work, represents a yearly 
mean mid-latitude temperature profile and does not include 
the seasonal or short-term variations in temperature nor does 
it describe well the arctic or tropical thermal structure. 
Listed in Table 19 is a comparison of the U.S. Standard 

V 

Atmosphere 1976 density profile to the mid-latitude summer and 
winter, the arctic summer and winter, and the tropical density 
profiles obtained from the. U.S. Standard Atmosphere 1962 (19). 

Although these .density variations are not entirely attri- 
butable to variations in the atmospheric temperature profile, 
the numbers in Tables 18 and 19 suggest that the error in the 
estimated optical airmass could be large if the temperature 
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profile, including the variation along the ray path, is not 
properly modeled. 

Inherent variations of the temperature profile both along 
the ray path and seasonally make the accurate retrieval of 
mixing-ratio profiles difficult. One method of avoiding this 
difficulty is to retrieve a temperature profile from the 
spectral data to be analyzed. Temperature dependent line 
intensities and the known mixing-ratio profile of gases such as 
CO 2 can be used to extract the thermal structure of the atmos- 
phere by analysis of the CO 2 lines in the spectra. It must be 
noted that the temperature profile thus obtained will be some 
mean profile over the entire path. Since, as discussed in 
section 2.4, the retrieved mixing-ratio profile is also a mean 
over the path, the mean temperature profile obtained from the 
experimental spectral data is preferable to the use of standard 
atmospheres . 

The previous discussion implies that these errors can be 
reduced if the temperature profile is determined internally 
(by using the same spectral data that is used to determine the 
mixing-ratio profile) . The relative advantages of internal 
calibration over the use of externally determined physical 
parameters is discussed later in this chapter. 

Uncertainty in the temperature profile 

The systematic errors due to the inherent error in the 
temperature profile can be estimated by a simple derivation. 
Assume that the atmosphere is isothermal and the pressure at 
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altitude A is given by 


P (A) = P (o) exp (-A/H g ) 


( 15 ) 


where H q , the pressure scale height is 


u _ R*T 
Hs “ Mi" * 


(16) 


From Eq. (15) at altitude A a fractional change in pressure 
AP occurs in an altitude change AA given by 


AP AA 
P ~H C 


(17) 


Equation (15) can be written in terms of the temperature by 
combining it with Eq. (16) yielding 


P (A) = P (o) exp (-AMg/R*T) 


:is) 


From Eq. (18) the change in temperature AT due to a change in 
the pressure AP is given as 


AP _ AMgAT 
P R*T 2 


(19) 


The ideal gas law gives the density at altitude A as 


p (A) = MP(A)/R*T . 


( 20 ) 
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A change in density Ap(.A) results in changes in the pressure 
and temperature given by 


Ap(A) 


MAP (A) 
R*T 


MP (A) AT 
R*T 2 


( 21 ) 


or 


Ap (A) _ AP (A) AT 
p (A) P (A) T 


Combining Eqs . (17) and (22) yields 


Ap (A) _ AAT AT 
p(A) H S T " T ' 


or 


Ap (A) 
P (A) 


(A - H„) AT 

O 

H^T 


(24) 


Since, from Eq. (10) , the fractional change in density equals 
the fractional change in the calculated optical airmass 
Eq. (24) becomes 


Aa (A -H S )AT 
a H S T 


(25) 


Thus, for an isothermal atmosphere uncertain by amount AT, 
the error in the estimated optical airmass will be zero at 
an altitude equal to the scale height. 
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The error due to temperature error in the U.S. Standard 
Atmosphere was investigated by comparing optical airmasses, 
corresponding to several tangent heights, calculated using the 
standard values of Table 15 to those obtained using identical 
values with the U.S. Standard Atmosphere profile incremented 
at all altitudes by ±5 K. Some of the results from this com- 
parison are listed in Table 20 together with the errors pre- 
dicted by Eq. (25) . 

As expected a change in sign of the error occurs near the 
scale height. The errors obtained for AT = ±5 K agree well with 
those predicted by Eq. (25) . The small differences are 
attributed to the non-isothermal structure of the atmosphere 
used. 

3.2.5 Type 3_ - observer altitude 

As shown in Figure 46 an error in the altitude AA Q of the 
observing instrument results in an error in altitude along the 
entire path of approximately AA 0 sinZ 0 or 

AA=AA 0 sinZ 0 . (26) 

Previously it was noted that the fractional decrease in 
pressure for an increase in altitude is given by Eq. (17) . 
Equivalently an increase in altitude AA will result in a 
decrease in pressure AP given by 
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Table 20: Error in computed optical airmass for systematic temperature error of ±5 

at all altitudes. 
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Figure 46: The effect of observer altitude error on the 

tangent height. 
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AP 

P 


Aa sinZ 
o o 


( 27 ) 


H s 

The optical airmass encountered along a ray path .increases 
proportionately if the pressure increases (see Eq.. (13)). 
Therefore 


Aa 

a 


AA sinZ 
o o 

H S 


(28) 


Thus for a particular apparent angle and an isothermal atmos- 
phere the fractional error in the estimated optical airmass 
traversed by a ray is proportional to the error in the obser- 
ver altitude. 

This relation was examined by comparing the optical air- 
mass calculated from the standard values of Table 15 to 
airmasses calculated using identical values but with different 
observer altitudes ranging from 499.5 km to 500.5 km. In 
these cases the maximum error in the airmass, from Eq. (28), 
should be 6.6%. Results obtained for three tangent heights 
are shown in Figure 47. The larger than expected error for 
the 21.96 km tangent height is probably due to departures of 
the standard atmosphere temperature profile from isothermal 
near the tangent height of this ray. The smaller error 
obtained for +AA 0 as compared to -AA 0 is believed to be due 
to refractive effects not considered in the derivation of 
Eq. (28). 
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AAo(KM) 

Figure 47: The error in the calculated optical airmass due 

to error in the observer altitude for tangent 
heights of 3.04, 21.96, and 43.77 km. 
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Thus for an observer altitude of 500 km, errors of 0.5 km 
in the altitude of observation result in errors of 6.6% in 
the calculated optical airmass. In order to maintain 1% 
error in the estimated optical airmass this altitude must be 
known to within 75 m, suggesting that high precision is 
necessary in determining the satellite height. 

3.2.6 Type 3_ - apparent zenith angle 

In practice the apparent zenith angle used in the slant 
path calculation is not usually measured directly. This angle 
can be determined from the geometry of the orbit by recording 
the time of observation and the portion of the solar disc 
viewed. If refraction is included the apparent zenith angle 
can then be calculated. Estimates of the effects of time 
error on the estimation of the optical airmass can be made 
from consideration of the errors in the apparent zenith 
angle. 

If there is an error AZ 0 in the measurement of the 
apparent zenith angle the true tangent height of the ray is 
shifted by AA^ (see Figure 48) . If AZ 0 is small then 

AA t --LAZ q (29) 

where L, the path length between the tangent height and the 
observer, is approximately 2500 km for a 500 km observer 
altitude. Combining Eqs. (13), (17), and (29) gives 
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Aa 

a 


L 


(30) 


AZ 0 

hT * 

Thus due to the long path lengths for an observer alti- 
tude of 500 km, an error of about 0.2 milliradians results 
in an error of about 7% in the airmass calculation. The 
errors due to the inherent error of the apparent zenith 
angle are large for any satellite experiment. 

The effects of apparent zenith angle errors were 
examined by comparing the optical airmass calculated for the 
standard values of Table 15 to those obtained using identi- 
cal values and an incremented apparent zenith angle. Some 
results, shown in Figure 49, were obtained for three dif- 
ferent tangent heights. The error predicted from Eq. (30) 
for an apparent zenith angle error of 0.01° is 6.2% for all 
three tangent heights. The variations seen in Figure 49 
between the three tangent heights , which were not predicted 
by Eq. (30) , are probably due to the temperature variation 
with altitude of the standard atmosphere. The asymmetry, 
about AZ o = 0, is believed to be due to the effects of refrac 
tion not included in this derivation. 

Smearing due to observer motion and the field of view 

For the case of a satellite-borne instrument, which 
moves in orbit during the time necessary to complete one 
occulta tion, the optical airmass encountered by rays 
reaching the instrument at the end of the collection period 
will be substantially different than that encountered by 
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Figure 49: The error in the calculated optical airmass due 

to error in the apparent zenith angle for tangent 
heights of 3.04, 21.96, and 43.77 km. 
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rays at the beginning of this period. Park (14) has shown 
that for the collection of an interferograxn this results in 
the information about the spectral line strengths and the 
spectral line half-widths being obtained from rays which 
traverse different portions of the atmosphere. In addition, 
because of the finite size of the solar disc viewed by the 
instrument, rays which reach the instrument at identical 
apparent zenith angles may have different tangent heights 
(15) . The effects of smearing of the spectral information 
due to motion of the observer and the finite field of view 
can be approximated by systematic errors in the apparent 
zenith angle. 

This can be demonstrated by use of Kepler's laws. We 
have assumed, in order to obtain an estimate of the smearing 
due to satellite motion, that the orbit of the satellite is 
circular, coplanar with the sun, a height A above a spherical 
earth of radius R e . If the source of radiation is assumed, 
for now, to be a point source at infinity, the change in 
apparent zenith angle AZ 0 of the incoming radiation is equal 
to the change in polar angle A0 of the satellite - , as shown 
in Figure 50. Thus AZ Q is given by 

AZ q = 2tt^ , (31) 

i-e 

where At is the time necessary for the satellite to move 
angular distance A0 and t e is the period of the satellite 
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Figure 50: The orbital geometry of a moving satellite 

observer. The apparent zenith angle of observa- 
tion changes by AZ 0 for a change A0 in the orbital 
position . 
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orbit. The period of the orbit is given by Kepler's laws. 
Therefore the rate of change of the apparent zenith angle 
is 


AZ q 

“St" 


/ GM e 

(R e + A 0 ) 3/2 


(32) 


where M e is the mass of the earth. 

The rate of change of the optical airmass encountered 
by rays reaching the instrument during an occultation was 
estimated by calculating the optical airmasses for apparent 
zenith angles determined from Eq. (32) for a constant time 
spacing At" = 0.1 sec. The optical airmass was computed for 
the standard values of Table 15. The fractional rate of 
change of the optical airmass was defined as 

a y ~ 3-n 

„ _ O (t) o(t + At") 

F — (33) 

Z o(t) 

where the units of F are sec/ The values of F(Z Q ) , shown 
in Figure 51, were obtained for the U.S. Standard Atmosphere. 
The rate of change of the airmass depends on the tangent 
height but is not less than 30%/sec. 

The structure of the F(Z 0 ) curve is due to the tempera- 
ture variation of the standard atmosphere, shown in Figure 40. 
The dependence of F(Z Q ) on the temperature profile was 
demonstrated by repeating the computations of Eqs. (32) and 
(33) using instead a temperature profile which decreased 
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Figure 51: The fractional time rate of change of the optical 

airxnass for a moving satellite observer and the 
standard atmospheric values as a function of 
apparent zenith angle. The tangent height range 
is 80 km to ground level. 
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linearly with altitude. The values of F(.Z 0 ), shown in 
Figure 52, were obtained for this temperature profile. The 
F(Z 0 ) determination was also made for the case of the U.S. 
Standard Atmosphere temperature profile but no atmospheric 
refraction. These values, shown in Figure 53, demonstrate 
that the shape of the F(Z 0 ) curve depends only on the atmos- 
perhic temperature profile. This suggests that, in order to 
estimate this rate of change and therefore the effect of this 
smearing on the spectral data, the temperature profile must 
be accurately known. 

In addition, for satellite occultation experiments such 
as ATMOS, the field of view of the instrument is 2 milliradians 
and thus, from Eq. (30), the variation in the optical airmass 
traversed by arriving rays can be as great as 70%. 

The effects of smearing due to the motion of the obser- 
ver and the finite field of view are larger than would be 
expected from errors in the apparent zenith angle alone, and 
these effects are not easily eliminated from the spectral 
data. However, the smearing effects on individual lines in 
a single spectrum are identical (for weak lines). Thus it 
is possible by comparison of spectral lines resulting from 
an absorber with a known mixing-ratio profile to spectral 
lines resulting from an absorber with unknown distribution 
to eliminate these weighting effects in the mixing-ratio 
profile determination of the unknown absorber. This method 
of internal calibration is discussed fully later in this 
chapter. 


166 



F(Z 0 ) 



n oc i 1 1 1 1 i 1 1 1 1 1 

no in 112 

APPARENT ZENITH ANGLE (DEGREES) 


Figure 52: The fractional time rate of change of the optical 

airmass for a moving satellite observer and a 
temperature profile which increases linearly with 
altitude as a function of apparent zenith angle. 
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RATE OF CHANGE OF AIR MASS 



Figure 53 : The fractional rate of change of the optical 

airmass, for a moving satellite observer and 
assuming no atmospheric refraction, as a function 
of apparent zenith angle. 
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3 .2.7 Type 3 - latitude 

The shape of the earth is a nearly spherical ellipsoid. 
Since the long atmospheric path lengths for a 500 km obser- 
ver altitude may result in latitude changes between the 
observer and the tangent height as large as 10° , it is impor- 
tant to know the magnitude of the errors resulting from changes 
in latitude along the ray paths. 

The change in optical airmass due to changes in the 
earth's latitude or equivalently the earth's radius can be 
estimated by starting with Fermat's principle (Eq. (A.l) of 
Appendix A) 

n 0 r 0 sinZ 0 = n-pr-p (34) 


where 


r T = A T + R e ‘ (35) 

If refraction is neglected Eq. (34) becomes 

(A 0 + R e ) sin Z Q = A T + R e . (3 6) 

A change in radius ARe results in a change in tangent height 
A Am given, from Eq. (36) , by 

AA<p = (sin Z 0 -l)AR e . (37) 
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Combining Eqs. (12) , (17) , and (37) yields 


Aa 

a 


(1 - sin Z Q ) AR e 

Hi 


( 38 ) 


The radius of curvature of the earth is given as a func- 
tion of latitude by (20) 


„ ... )c 4 cos^d) + b 4 sin^<j> 

RM<M 5 2 o 

c cos^tp + b sin^ 


(39) 


where 


c = 6378.160 km 


(40) 


and 


b = 6356.775 km . (41) 

Thus the fractional error in the airmass can be expressed as 
a function of latitude by using Eq. (38). The change in 
optical airmass for a latitude change from 45° to 55° is 
found to be 3.3%. The maximum error in using a latitude of 
45° rather than the true latitude is found to be 10%. Though 
not large in magnitude for reasonable latitude errors (±5°) 
higher precision retrievals may be obtainable if the effect 
of the changing earth's radius is modeled. 
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The validity of Eq . (38) was tested by varying the 
earth l s latitude with respect to 45°. The error obtained 
in the optical airmass for the standard value of Table 15 
is shown in Figure 54 for rays with four tangent heights. 

These errors are much larger than those predicted by Eq. (38), 
although the shapes of the curves agree with that predicted 
by Eq. (38).. Figure 55 shows ARq as a function of latitude, 
which exhibits the same dependence as the error curves of 
Figure 54. Similar error curves, but with a much smaller 
amplitude, shown in Figure 56, were obtained for an observer 
altitude of 50 km. The large dependence of the amplitude on 
the observer height suggests that the amplitude of the error 
curves is incorrect although the functional form is consistent. 
It is believed that the large amplitude for observers above 
50 km is due to an error in the tangent height calculation 
algorithm used in the slant path calculation (8) . This error 
is not expected to alter any of the o'ther results obtained 
in this work. 

3.2.8 Summary 

In this section we have estimated the errors introduced 
into the determination of the optical airmass traversed by 
a ray passing through the atmosphere due to improper layer 
thickness in the calculation, neglect of the effects of 
refraction, errors in the sea-level pressure, the atmospheric 
temperature profile, the observer altitude, the apparent 
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LATITUDE OF OBSERVER (DEGREES) 


Figure 54 : The error in the calculated optical airmass due 

to error in the latitude with respect to 45° for 
tangent heights of 2.63, 5.86, 15.01, and 51.19 km 
and an observer altitude of 500 km. 






LATITUDE (DEGREES) 


Figure 56: The error in the calculated optical airmass due to 

error in the latitude with respect to 45° for 
tangent heights of 25.5, 41.2, and 49.0 km and an 
observer at 50 km. 
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zenith angle of arrival of the ray at the observer, and the 
latitude. By far the largest systematic error in the opti- 
cal airmass calculation was caused by error in the apparent 
zenith angle. 

A much more interesting problem is that of the effects 
of these systematic errors on the retrieval of mixing-ratio 
profiles from absorber amount data. This problem is considered 
in the next section. 

3.3 THE EFFECT ON THE MIXING-RATIO RETRIEVAL 
3.3.1 A general formalism of systematic errors 

The amount of absorber a^ encountered by radiation 
traversing the atmosphere is given by. 

N 

a m = I m^pjj65,^ . (42) 

j=l J 

It is assumed in the following analyses that the absorber 
amount can be exactly determined from analysis of experimen- 
tal data (i.e., equivalent width analysis) and that Aa m = 0 . 

This implicitly assumes that the absorber amount data 
considered here contain no experimental noise. The effects 
of noise, described in the previous chapter, must be added to 
the effects due to systematic errors examined in this section 
in order to determine the total error in the retrieved 
mixing- ratio profile. 
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Since Aa m = 0, the variation of Eq. (42) due to variable 
x becomes 

N Ap 

0 = . Z 1 (P j <S£ j 13T Ax + m j 6ii j “5x Ax + m j p j 1 Ax) 

(43) 

If Ax represents a change in the air density Ap and if the 
effects of refraction are ignored then the A6£j/Ax terms 

are identically zero since the ray path is a straight line 
from source to observer. In this case, since the retrieved 
mixing-ratio values are assumed to be linearly indepednent, 
Eq. (43) reduces to 


p j im j + m j Ap j = 0 


or 




But from Eq. (6) 


Ap j _ Ap (o ) 
Pj P (o) ' 


so 



Ap (o ) 
P (o) 


(44) 


(45) 


(46) 


( 47 ) 


176 



which combined with Eq. (10) yields 



Consequently, the error in the retrieved mixing-ratio value 
at any altitude is given by Eqs . (13), (25), (28), (30), and 

(38) . 

The effects of systematic errors in the parameters in 
these relations on retrieved mixing-ratio profiles have been 
investigated by analyzing absorber amount data calculated for 
the standard values of Table 15 by the methods, described in 
the previous chapter assuming various systematic errors in 
the calculation of the a-^j matrix elements of Eq. (2.70) . 

The 125 absorber amounts were calculated for the constant 
mixing-ratio profile of Figure 11 and an apparent zenith 
angle spacing of 0.0098903°. In these retrievals no simu- 
lated noise was added to these data values. The effects of 
the systematic errors were investigated by comparing mixing- 
ratio values for Nr = 50 to the true mixing-ratio profile. 

3.3.2 Type 1 - refraction 

If refraction is ignored in the calculation of the 
retrieval matrix substantial errors occur in the retrieved 
mixing-ratio values at low altitudes. The error in the 
retrieved mixing-ratio values as a function of altitude due 
to the neglect of refraction is shown in Figure 57. Below 
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25 km the error increases rapidly near the earth's surface. 
Above 25 km the error in the retrieved mixing-ratio values 
is less than 1%. Thus refraction can be ignored above 25 km 
but must be properly modeled below 25 km in order to avoid 
systematic errors in the retrieval. The feature at 10 km, 
shown in Figure 57, is believed to be due to the shape of the 
temperature profile at the tropopause and the thickness of 
the layers used in the calculation. 

3.3.3 Type 2_ - sea-level pressure 

The error in the jth retrieved mixing-ratio value is 
found by combining Eqs . (13) and (48) yielding 

= Ap (o) (49) 

m j P(o) * 

Thus the error in the retrieval should be independent of 
altitude. For an error of 5%' in the sea-level pressure the 
error in the mixing-ratio retrieval should be 5% at all 
altitudes . 

Errors in the retrieved mixing-ratio profile obtained 
for sea-level pressure errors of ±5% are shown in Figure 58. 

As expected the error is about ±5% at all altitudes. The 
departure from. this, at low altitudes is due to the effects of 
refraction which have not been included in the derivation of 
Eq. (49) . The feature in the error curves at 10 km is believed 
to be due to the shape of the temperature profile at the 
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Figure 58: The error in the retrieved constant mixing-ratio 

profile due to errors in the sea-level pressure 
of ±5%. 





tropopause and the layer thickness used for the calculation. 

3.3.4 Type 2 _ - temperature profile 

The error in the jth retrieved mixing-ratio value for an 
error in an isothermal temperature profile of AT at all alti- 
tudes is found as a function of altitude by combining Eqs. (25) 
and (48) giving 


Anij 

m ■ ~ 


(A - He) AT 


H S T 


(50) 


Thus the error in the mixing-ratio value at an altitude equal 
to the scale height is zero independent of the temperature 
profile error. At 50 km the error in the retrieval is ±11% for 
a systematic temperature error of ±5 K at all altitudes. The 
problem of local temperature variations has not been considered. 

The results of retrievals done for systematic temperature 
errors of ±5 K are shown in Figure 59. As expected the error 
in the mixing-ratio retrieval is zero near an altitude equal to 
the scale height and increases in magnitude with altitude above 
this point. The error at 50 km is larger than that predicted 
by Eq. (50) because this equation was derived assuming an 
isothermal atmosphere. 

The altitude dependence of the mixing-ratio value error 
caused by temperature profile errors is unique to this syste- 
matic error (of those discussed in this work) . Thus the 
inherent error in the retrieved mixing-ratio profile due to 
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Figure 59: The error in the retrieved constant mixing-ratio 

profile due to errors in the temperature profile 
of ±5 K at all altitudes. 





temperature errors cannot be reduced by adjusting the other 
parameter values. Therefore, the temperature profile should 
be known accurately and, if necessary derived from the same 
spectral data used for the mixing-ratio profile retrieval (21) . 
Although not discussed in this work, it is planned to describe 
a method to retrieve temperature profiles later. 

3.3.5 Type 3_ - observer altitude 

The error in the retrieved mixing-ratio value due to 
error in the observer altitude is found by combining Eqs . 

(28) and (48) giving 

Am. AApsin Z Q (51) 

mj 

This error, like that obtained for sea-level pressure errors, 
should be independent of altitude. For an observer altitude 
error of 1 km the error predicted from Eq. (51) is ~"14% at all 
altitudes . 

The mixing-ratio retrieval value errors, shown in Figure 
60, were obtained for an error of ±1 km in observer height 
from analysis of the absorber amounts, calculated from the 
standard values. The mean value of the error is approximately 
14% as predicted. However, the features shown in Figure 60 
are believed to be due to the temperature variation of the • 
standard profile and the layer thickness used in the calcula- 
tion. 
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Figure 60 : The error in the retrieved constant mixing-ratio 

profile due to errors in the observer altitude 
of ±1 km. 
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3.3.6. Type 3_ - apparent zenith angle 

Similarly the error, due to errors in the apparent 
zenith angle, in the jth layer is given by 




( 52 ) 


An apparent zenith angle systematic error of 0.01° results in 
an error in the retrieved mixing-ratio value of 6% at all 
altitudes . 

The results of retrievals done for apparent zenith angle 
errors of ±0.01° are shown in Figure 61. As expected the 
error is nearly constant with altitude and has a mean value of 
approximately 6%. The structure in the error curves is almost 
identical to that obtained for the observer altitude errors. 
The similarity of Figure 60 and 61 suggests that systematic 
errors in the apparent zenith angle are indistinguishable from 
systematic errors in the observer altitude. 


3.3.7 Type 3_ - latitude 

The mixing-ratio error due to latitude error is given by 
Am • (1 - sin Z Q ) AR e (<j>) 



For a latitude error of ±5° with respect to the standard value 
of 45° the error predicted by Eq. (53) is 1.8%. This error 
like that found for sea-level pressure, observer altitude, and 
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Figure 61: The error in the retrieved constant mixing-ratio 

profile due to errors in the apparent zenith angle 
of ±0.01°. 




apparent zenith angle errors should be independent of altitude' 
in the absorbing atmosphere. 

The results obtained from the analysis of the standard 
absorber amount data, shown in Figure 62, gives an error three 
times larger than that predicted by Eq. (53). This is believed 
to be due to the algorithmic error in the slant-path calcula- 
tion mentioned in the previous section. The shape of the 
error curves, however, are identical to that obtained for the 
observer altitude and apparent zenith angle errors. 

3.3.8 Summary 

By analysis of synthetic absorber amount data with no 
added noise, assuming various systematic errors in the de- 
scribing parameters, we have determined the error in retrieved 
mixing-ratio profiles as a function of altitude caused by 
these systematic errors . 

It was found that refraction becomes significant in the 
retrieved profile below 25 km and its importance increases 
rapidly with decreasing altitude. In order to retrieve 
accurate mixing-ratio profiles below 25 km the effects of 
atmospheric refraction must be modeled properly. 

The error in the retrieved mixing-ratio profiles obtained 
for systematic errors in the temperature profile increased 
linearly with altitude and was zero at an altitude approximately 
equal to the scale height. The functional form of this syste- 
matic error was different from that obtained for all the others 
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Figure 62: The error in the retrieved constant mixing-ratio 

profile due to errors in the latitude of ±5°. 
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suggesting that this error can only be reduced by a more precise 
determination of the temperature profile from analysis of the 
same spectral data used to retrieve the mixing-ratio profiles. 

The errors due to systematic errors in the sea-level pres- 
sure, the apparent zenith angle, the observer altitude, and 
the latitude are essentially identical and independent of 
altitude. Consequently, the effects of an error in any of 
these parameters on the retrieval of a mixing-ratio profile 
are indistinguishable from the effects of errors in the other 
parameters. This implies that a systematic error in the 
observer altitude could be compensated by an appropriate adjust- 
ment of the other parameters. The result would be the retrieval 
of a mixing-ratio profile with no apparent systematic errors 
(provided that refraction is modeled and the temperature pro- 
file is precisely known) . A method of internal calibration of 
spectral data which makes use of this technique will be 
described later in this chapter. 

3.4 ESTIMATE OF TOTAL SYSTEMATIC ERROR 

Although individually the contribution of each of the 
systematic errors discussed in the previous section may be small 
compared to the uncertainty caused by the presence of experi- 
mental noise, in combination these errors can significantly 
alter the results obtained in a mixing-ratio profile retrieval. 
For the design of occultation experiments it is necessary to 
know the accuracy required for the determination of the physical 
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parameters in order to maintain a desired error in the mixing- 
ratio profile retrieval. Conversely, it may also be necessary 
to determine the total error due to possible errors in these 
physical parameters. Specifically, for the evaluation of 
experimental results, the total error in the retrieved mixing- 
ratio values due to the presence of systematic error can be 
computed. 

In this section the experimental design requirements for 
an error in trie retrieved mixing-ratio values of 1% for each 
of the systematic errors (excluding refraction) discussed in 
the previous section is investigated. This represents a total 
systematic error of about 5% in the retrieved mixing-ratio 
values. In addition, the altitudes below which each of the 
errors is significant compared to the uncertainty in the 
mixing-ratio value due to experimental noise is computed. The 
calculation of this altitude allows the determination of 
those portions of the retrieved profile significantly affected 
by the systematic errors. 

In order to determine the magnitudes of the systematic 
errors necessary to give a total error of 5% Eqs. (49) -(53) 
were solved for each error by substituting in Amj/mj=0.01. 

The parameter errors which result in a mixing-ratio value 
error of 1% are listed in Table 21. These systematic errors 
have been listed in order of the accuracy necessary to obtain 
an error of 1%. For example, an apparent zenith angle error 
of 0.002° is a change of only 0.002% with respect to the 
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Table 21: Systematic errors in physical parameters which give 

a 1% error in the retrieved mixing-ratio profiles. 


Physical 

parameter 

Dependence of error 
in mixing-ratio on 
height 

Systematic error 
which gives a 1% 
error in mixing- 
ratio profile 

refraction 

important below 25 km 

at all altitudes 
below 25 km if 
neglected 

apparent zenith 
angle 

independent of height 

0.002° 

satellite height 

independent of height 

75 m 

temperature profile 

increases with height 

0.1 K at 25 km 

sea-level pressure 

independent of height 

10 mb 

latitude 

independent of height 

2.5° 
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standard values of Table 15, while a latitude error of 2.5 
represents a change in latitude of about 5.6% with respect to 
the standard value. Thus the error in the apparent zenith 
angle was judged the more critical, and so on. 

The errors in the experimental parameters listed in 
Table 21 indicate the precision necessary in an experiment to 
yield a total systematic error of 5% in the retrieved mixing- 
ratio values. It is not likely, due to smearing caused by the 
motion of the observer and the finite size of the solar disc 
viewed by the instrument, that an overall precision as high as 
this is possible for satellite experiments. The magnitudes 
assumed in the previous section for the systematic errors are 
more probable. 

The magnitudes of the systematic errors used in the pre- 
vious section are listed in Table 22. These errors may result 
in a total error in the retrieved mixing-ratio values of up 
to 30%. A measure of the significance of these systematic 
errors is the altitude range over which these errors are larger 
than the uncertainty in the retrieved mixing-ratio values due 
to the uncertainty of experimental noise. These altitudes were 
calculated, and are listed in Table 22, for each of the syste- 
matic errors and the conditions: 125 data values with an SNR 

of 100 and the retrieved number of mixing-ratio values equal 
to 50 and 10. Thus for a mixing-ratio value spacing of 10 km 
the systematic errors are individually larger than the uncer- 
tainty due to noise below 24 km. In this case more than half 
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Table 22: The altitude at which the error in the mixing-ratio 

profile due to each systematic error is less than 
the uncertainty due to noise for N 0 = 125, 

SNR = 100, and Nr values of 50 and 100. 


Physical 

parameter 

Systematic 

error 

Altitude at which 
systematic error is 
less than uncertainty 
due to noise 



Nr = 

50 

N r = 100 

apparent zenith 
angle 

0.01° 

6 . 6 

km 

17 km 

satellite height 

1 km 

12.0 

km 

23 km 

temperature profile 

5 K 

0.0 

km 

4 km 

sea- level pressure 

5% 

5.4 

km 

16 km 

latitude 

5° 

0.0 

km 

9 km 
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the altitude range of the retrieval will contain mixing- 
ratio value errors primarily due to systematic errors. Thus, 
in order to obtain meaningful profile retrievals it is neces- 
sary that these errors be reduced. A method for reducing the 
effect of these systematic errors is discussed in the next 
section. 

The relations derived in this chapter for the retrieved 
mixing-ratio value errors due to systematic errors are valid 
for all observer altitudes. The errors due to observer alti- 
tude, sea-level pressure, and temperature profile errors are 
essentially the same for all observer altitudes. For the case 
of a balloon observer (i.e., an observer altitude of about 30 km) 
the error in the retrieved mixing-ratio values caused by an 
error in the apparent zenith angle (for apparent zenith angles 
greater than 90°) of 0.01° is about 2%, and the error due to an 
error in the latitude of 5° is 0.1%. Thus the systematic error 
due to these parameters is much less for a balloon observer. 
However, as shown in Table 23, typical errors in the physical 
parameters give a total error of about 14% in the retrieved 
mixing-ratio values independent of altitude. This suggests that 
many of the reported mixing-ratio profiles retrieved from the 
analysis of data obtained from balloon-borne instruments may 
have been overly optimistic in their reported uncertainties. 
Realistic errors in the physical parameters are necessary if 
the accuracy of retrieved profiles is to be gauged. 
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Table 23: The error in retrieved mixing-ratio profiles due 

to typical systematic errors for an observer 
altitude of 30 km. 


Physical 

parameters 

Systematic 

error 

Error in mixing-ratio 
profile 

apparent zenith 
angle 

0.010 

2.00% 

satellite height 

0.5 km 

7.00% 

temperature profile 

5 K 

0.25% 

sea-level pressure 

' 5% 

5.00% 

latitude 

1° 

0.02% 
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For the case of a ground-level observer the systematic 
errors due to the apparent zenith angle, latitude, and observer 
altitude errors are eliminated. The remaining errors due to 
errors in the sea-level pressure and the temperature profile 
(including errors due to "horizontal" variations of the sea- 
level pressure and the temperature profile) should be small 
compared to the error in the retrieved mixing-ratio profile 
caused by noise. Thus it is not expected that these systematic 
errors will have much affect on the retrieval of mixing-ratio 
profiles from ground-level observations unless the SNR of the 
data is large. 

3.5 INTERNAL CALIBRATION OF OCCULTATION SPECTRA 
3.5.1 Internal calibration method 

It was demonstrated in section 3.3 that when the informa- 
tion about the atmospheric paths is obtained from external 
sources several errors affect the retrieved profiles in nearly 
identical ways. This occurs in part because the spectral 
information concerning the atmosphere, carried by a ray with 
a fixed tangent height originating outside the atmosphere and 
collected outside the atmosphere, is not influenced by the 
positions of the source or the detector and hence by the 
particular sets of parameters which describe these positions. 

As shown in Figure 63 the ray path through the atmosphere is 
identical for an apparent zenith angle of 110.807° together 
with an observer altitude of 500 km to the path of a ray with 
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apparent zenith angle of 110.587° together with an observer 
atltitude of 490 km. If the tangent height of the observed 
radiation can be determined the particular values assumed for 
the observer parameters do not affect the spectral analysis as 
long as the tangent height corresponds to that obtained from 
the data. Thus it is possible to choose pairs of values for 
the apparent zenith angle and the observer altitude such that 
the ray path is consistent with that obtained from the data. 

This is the essence of the internal calibration method. 

It has also been noted that the quality of the spectra 
is degraded by the field of view of the instrument, satellite 
motion (14) , atmospheric refraction, and temperature and 
density profile variations along the atmospheric path (15) and 
that it is unlikely that all of these effects can be satis- 
factorily modeled in the analysis. , 

Internal calibration removes many of these errors but also 
introduces new ones. In this method information about the 
atmospheric path is obtained from absorption lines in the 
spectra. Lines of CC >2 are often measured because its vertical 
and latitudinal mixing-ratio is nearly constant and the para- 
meters of many CC >2 lines are accurately known. We consider, 
initially, the information obtained by measuring the equivalent 
widths of a sing'le CO 2 line, with a temperature-independent line 
intensity, in a single spectrum. This allows one aspect of the 
path, such as the total amount of CO 2 (which gives the tangent 
height of the ray) , to be obtained. The additional information 
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obtained by measuring the equivalent width of several lines 
in different regions of the curve of growth or by measuring 
lines whose intensities are temperature dependent is considered 
later in this section. 

If W Q and W]_ are the equivalent widths, in a single spec- 
trum, of weak lines of CO 2 and a trace gas, respectively, and 
if the corresponding line intensities S Q and Sj_ are temperature 
independent, then from Eq. (2.9.) the amount of CO 2 along the 
path is 


ao — W 0 /S 0 

and the amount of the other absorber is 


(54) 


a l = w l/ s l • 


(55) 


If the CO 2 mixing-ratio m 0 is constant with height and 
along the ray path traversed then the total optical airmass 
traversed is 


a 


w 


^o m o 


(56) 


and the mean mixing-ratio of the absorber along the path is 
give by 


mi 


a l 


a ’ 


(57) 
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These simple results are obtained because the incremental 
equivalent widths are additive. In this approach it is assumed 
that both equivalent widths have been subjected to the same 
averaging due to the smearing caused by the field of view of 
the instrument, the satellite motion, and the temperature and 
density profile variations. If the spectra are obtained by 
a Fourier Transform Spectrometer information concerning both 
of the lines W Q and W]_ is collected simultaneously. Provided 
the lines are both weak, both of these equivalent widths are 
subjected to the same weighting processes which are then auto- 
matically ratioed out of the final results. Thus the values 
or the total airmass traversed and the mean mixing-ratio of the 
unknown absorber given by Eqs. (55) and (57) do not contain the 
implicit errors associated with the uncorrectable smearing 
effects of the external calibration methods. 

The elimination of the smearing effects represents the 
principle advantage of the internal calibration method. These 
effects are extremely difficult to characterize analytically, 
yet by internal calibration of the spectral data this difficulty 
is eliminated. The problem then reduces to the idealized 
analysis described in the previous chapter. 

The corresponding effective tangent height of the rays pro- 
ducing the spectrum can be found from comparison of the total 
airmass traversed to the schematic optical airmass plot, shown 
in Figure 64, which could be obtained from the slant-path cal- 
culating program (7,8). The tangent height for an optical 
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airraass a is found from Figure 64 to be A T . 

If we consider N c spectra obtained for different paths 
through the atmosphere, the N 0 tangent height values could be 
obtained from Figure 64. A set of effective zenith angles 
together with an effective observer altitude can be chosen to 
give the correct tangent heights given from Figure 64. The 
set -of constituent absorber amounts (obtained from equivalent 
width analysis of a spectral line of the constituent) together 
with the corresponding effective apparent zenith angles and 
observer altitude constitute a data set identical to those 
considered previously in this work. However, the systematic 
errors due to the apparent zenith angle and the observer alti- 
tude errors as well as smearing due to satellite motion, field 
of view of the instrument, and temperature and density profile 
variations have been eliminated by the internal calibration 
method . 

In order to calculate the optical airmass curve shown in 
Figure 64 it is necessary to assume values for the temperature 
profile and the reference level pressure. If either of these 
is erroneous, then from Eqs. (13) and (25) the optical airmass 
curve of Figure 64 will be shifted from its true position. 

Thus the tangent height for each of the N 0 spectra will be 

improperly determined. However, the retrieved mixing-ratio . 

values will not be in error since the density of the incorrect 
density profile at this altitude is the same as the correct 
density at the correct tangent height. 


202 



Other systematic errors result if the line intensities 
of the CO 2 lines, the mixing-ratio value of CO 2 or the 
measured equivalent widths of the CO 2 lines are incorrect. 

The optical airmass determined using erroneous values for any 
of these parameters will also be in error, and consequently, 
the corresponding tangent height and the retrieved mixing- 
ratio values. Finally, although systematic errors in the 
constituent line intensities and the measured constituent 
equivalent widths do not affect the determination of the tan- 
gent heights of the rays, these errors result in errors in the 
constituent absorber amount and consequently, the retrieved 
constituent mixing-ratio values. 

These systematic errors are summarized in Table 24 together 
with the errors in the retrieved mixing-ratio values due to 
typical systematic errors. For the case of a 2% error in each 
of the spectral line intensities and measured equivalent widths 
the total error in the retrieved mixing-ratio values will be 
about 5% at all altitudes. No errors in these retrieved mixing- 
ratio vaules result from systematic errors in either the refer- 
ence level pressure or the temperature profile. However, these 
systematic errors as well as those due to the C0 2 spectral line 
parameters cause a shift in the altitudes of the retrieved 
mixing-ratio values due to the error in the tangent heights 
obtained from Figure 64. For the typical errors listed in 
Table 24 the total shift in the retrieved mixing-ratio profile 
should be about 0.4 km. 
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Table 24: The effects of systematic errors on mixing-ratio 

profiles retrieved using internal calibration. 


Physical 

parameter 

Systematic 

error 

Error in 
tangent height 

Error in 
mixing-ratio 
of trace gas 

surface 

pressure 

3% 

0.20 km 

none 

temperature 

profile 


depends on 
altitude 

none 

C0 2 intensity 
So 

2% 

0.14 km 

2% 

C0 2 equivalent 
width W Q 

2% 

0.14 km 

2% 

CO 2 mixing- 
ratio m Q 

2% 

0.14 km 

2% 

trace gas 
intensity 

2% 

none 

2% 

trace gas 
equivalent 
width W]_ 

2% 

none 

2% 
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An error of 0.4 km in the tangent height is much less than 
the expected vertical resolution of trace gas profiles expected 
to be obtained from typical occultation data. As shown in the 
second chapter, this resolution is typically several kilometers 
with uncertainties in the profiles of 10% or greater. Thus the 
method of internal calibration appears to be adequate for these 
types of measurements . 

The internal calibration method is useful in eliminating 
the systematic errors due to observer geometry parameter and 
atmospheric thermal structure parameter errors. This method 
also eliminates the smearing effects mentioned previously by 
comparison of constituent equivalent widths to those obtained 
from a gas with a known constant mixing-ratio profile such as 
CO 2 . The additional systematic errors which result from errors 
in the spectral line parameters are small enough to imply that 
internal calibration methods should be used to analyze spectra 
which contain suitable weak lines. 

3.5.2 Relative times of data collection 

The rate of change dZoi/dt of the zenith angles of the 
incoming rays and the corresponding rate of change dAip/dt of 
the tangent height of the rays reaching a satellite are nearly 
constant in a given occultation. Thus the scatter in the 
observed values of the optical airmasses obtained from Eq. (56) 
or the equivalent tangent height A^i with time allows the con- 
sistency of the data to be estimated. Systematic shifts between 
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the values obtained by internal calibration of the spectra and 
those calculated from the orbital parameters may be due, for 
example, to timing errors. Differences in the shapes of the 
curves at low tangent heights may be caused by improperly 
corrected refraction effects. However, the data in Figure 59 
indicate that small errors in the temperature profile do not 
cause significant errors, due to refraction, in the mixing-ratio 
profiles . 

Uncertainties in the CC >2 mixing-ratio or departure from 
local thermodynamic equilibrium may produce differences between 
the observed and calculated curves for small airmasses. 

However, good agreement is expected in the shapes of the curves 
for rays with tangent heights between about 20 and 50 km. 

3.5.3 Measurements of several lines 

The bands and lines of CO 2 have a wide range of intensities 
and it is expected that the equivalent widths of a number of 
lines can be measured in each spectrum. This enables the con- 
sistency of the measurements in a given spectrum to be estimated 
and allows those lines which are too weak to be measured 
accurately, which lie outside the linear region of growth of the 
equivalent width, or which are blended with other lines to be 
rejected. The effective SNR can be improved by measuring 
a number of lines. 

Because of the very large variation in the airmasses 
traversed by rays with tangent heights from 0 to above 50 km 
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it is necessary to measure CO2 lines with a corresponding range 
of temperature independent intensities in order to maintain high 
precision. In general, this will require measurements of CO2 
lines, with small lower state energies, in bands with different 
intensities . 

3 . 5.4 Lines not in the linear region 

It has so far been assumed that the equivalent widths of 
the CO2 calibration lines can be obtained from 

W(Z 0 i) = EAW-j_ j . ( 58 ) 

j 

If the CO 2 lines are not in the linear region of growth 
the incremental widths AWj_j depend on the local atmospheric 
temperature and density and they are no longer additive. The 
weighting of such lines over the field of view of the instrument 
and other smearing effects is different from that for weaker 
CO2 lines or weak lines of trace gases. 

In these cases one of the principle advantages of the 
internal calibration method, the ability to assume that the 
smearing effects are the same for all lines, is lost. It is 
therefore, important to choose CO2 calibration lines which are 
in or close to the linear region of the curve of growth.. 
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3.5.5 Measurements of temperature dependent lines 

Because the atmosphere is not isothermal the weighting 
of the equivalent widths of the lines over the field of view 
depends on the temperature dependence of the lines in addition 
to other effects previously considered. Since the strongest 
lines of most trace gases are not strongly temperature depen- 
dent, CC >2 lines with a similar temperature dependence should 
be chosen for calibration. 

3 . 6 SUMMARY 

In this chapter we have, investigated the effects of various 
systematic errors on the determination of the optical airmass 
traversed by a ray passing through the atmosphere and the 
mixing-ratio profiles retrieved from analysis of absorber amount 
data containing no noise. Expressions have been obtained for 
each of these resulting errors in the optical airmass and the 
retrieved mixing-ratio values as a function of altitude. 

It was found that refraction becomes significant below 
25 km and its importance increases rapidly with decreasing 
altitude. In order to retrieve accurate mixing-ratio profiles 
below 25 km the effects of atmospheric refraction must be 
modeled correctly. 

The. error in the calculated optical airmass and the 
retrieved mixing-ratio values due to systematic errors in the 
temperature profile of the entire atmosphere increased linearly 
with altitude and was zero at an altitude approximately equal 


208 



to the scale height. 

The errors due to systematic errors in the sea-level 
pressure, the apparent zenith angle, the observer altitude, 
and the latitude are essentially identical and independent 
of altitude. However, the errors, due to systematic errors 
in the latitude were found to be negligible for reasonable 
errors in the latitude. 

A method was presented by which the systematic errors due 
to temperature profile, sea-level pressure, apparent zenith 
angle, and observer altitude errors can be eliminated in the 
retrieval of mixing-ratio profiles from equivalent width data 
by compaxison of the equivalent width of a spectral line of the 
trace gas to the equivalent width of a CO 2 line in the same 
spectrum. If the C0 2 lines used for comparison are weak, this 
internal calibration method also eliminates the smearing effects 
due to the field of view of the instrument, satellite motion, 
and temperature and density profile variations along the ray 
paths. The systematic errors inherent to the internal cali- 
bration method should be less than those resulting from external 
calibration methods. 

In this chapter we have not considered the information 
contained in the shape of spectral lines. This problem is 
considered in the next chapter. Als o considered in the next 
chapter is the analysis of ATMOS spectra obtained from a recent 
ground-level test of the ATMOS instrument. As noted in this 
chapter the systematic errors caused by errors in the observer 
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geometry parameters are not significant for the analysis of 
spectra obtained from ground-level providing the SNR of the 
data is not large. Thus systematic errors due to errors in 
these parameters are not considered in the next chapter. 
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Chapter IV 

SPECTRAL CURVE-FITTING 
4.1 INTRODUCTION 

In the second chapter the problem of the retrieval of 
vertical mixing-ratio profiles from the analysis of absorber 
amount data was discussed. The information content of such 
data was estimated by determining the uncertainty and verti- 
cal resolution of the absorber mixing-ratio profile as 
a function of altitude. In the third chapter the effects of 
systematic errors on this profile were estimated by assuming 
incorrect values for the required geometric and atmospheric 
parameters. The data analyzed so far consisted of one piece 
of information, the total amount of absorber along the path 
(equivalent width data) , for each spectrum collected. In 
this chapter- the retrieval of mixing-ratio profiles by fitting 
the shape of spectral features is considered. 

The problem of spectral curve- fitting is treated in two 
parts. The first part is the retrieval of a constant mixing- 
ratio profile from a single spectral line for a ground-level 
observer. In the weak line approximation this reduces to 
a linear inversion problem to which the analysis techniques 
described in the previous two chapters can be applied. The 
uncertainty as a function of altitude was obtained for 
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various mixing-ratio value spacings. 

In the second part the results of the error analysis in 
the previous chapters is used to analyze spectral data 
obtained from a ground- level test of the ATMOS instrument. 
Non-linear least squares analysis (1) was used to retrieve 
mixing-ratio profiles for N 2 O and CH 4 from these data. 

These retrievals demonstrate that as expected there is 
little height discrimination in a single ground-level occupa- 
tion spectrum even if the spectral line shapes are modeled. 


4.2 THE INFORMATION CONTENT OF SPECTRAL LINE SHAPES 
4.2.1 General comments 

The data analyzed in the previous chapters were the sums 
of the absorber amounts in individual layers 

^ij = ra jPj^^ij * (1) 

Each of the a m — consisted of the desired mixing-ratio weighted 
by the local air density and path length. The ability 'to 
discriminate the vertical distribution of the absorber depended 
on the variations in the path segment lengths and on the 

angle of arrival of the rays. Since the weighting of the 
mixing-ratio by the air density is the same for all angles Z Q 
the density, as such, does not influence the derived vertical 
resolution directly. It does, however, affect the profile 
uncertainty. Indeed, the curves in Figures 16, 18, and 19 show 


214 




that the mixing-ratio uncertainty increases by, a factor of 
100 for a 30 km increase in altitude, both above and below the 
observer, in a manner which is directly proportional to the 
variation in the air density. 

Goldberg (2) has shown that, for rays with Z o <90° and no 
refraction, the path elements 

5£ = (1 + A/Re) <5A/(cos 2 Z 0 + 2A/R e + A 2 /R e 2 ) ** (2) 

where Rq is the radius of the earth and 6A is the layer thick- 
ness. As Re-*-* 

5£=secZ 0 6A . (3) 

In this plane-parallel approximation there is no change in the 
relative lengths of path segments through layers at different 
altitudes with incoming angle Z Q and hence no height discrimi- 
nation can be obtained. 

For altitudes below 50 km in the earth's atmosphere Eq. 02) 
can be approximated by 

61 - (cos 2 Z 0 + 2A/R e ) ~^6A . (4) 

Significant differences between Eqs. (3) and (4) occur as 
Z o -*-90° and the minimal resolution above the instrument, illus- 
trated in Figure 21, is due to small deviations from this 
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behavior... However, at these large angles, other effects such 
as atmospheric refraction, smearing caused by the portion of 
the solar disc viewed, and the time required for an observation 
must also be carefully considered if accurate profiles are to 
be retrieved. 

All the rays reaching an instrument, located above the 
surface, at angles greater than 90°, pass the earth at some 
minimum distance D. If refraction is neglected, the path length 
through a layer a distance A above this minimum is 

SSL-l (Re + D)/2 (A - D) ] ^6A -, (5) 

or, with respect to the minmum height, 

6 Z- VR^/2A') ^6A , (6) 

where R e ^ = R e + D and A" = (A - D) . 

The changes in the lengths of the path segments 6?,, and 
hence the weighting due to the geometry, are very large near 
the minimum height which, in turn, varies rapidly with Z Q . 

Some typical weighting curves have been described by 
Farmer (3) and Conrath (4) . The weighting curves for the 
analysis of absorber amount data in the second chapter decreased 
exponentially with altitude. It is these large variations 
which allow the high resolution, described by Eq. (2.83), to 
be obtained below the observer. 
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Since the geometry associated with occultations observed 
from the earth and the physical characteristics (especially 
the density profile) of the atmosphere are nearly constant, 
then the vertical resolution and uncertainty for a given pro- 
file obtained from equivalent width data with a fixed SNR and 
N 0 are completely determined. Any additional improvement must 
therefore come from the measurement of other spectral proper- 
ties . 

4.2.2 Vertical profiles from line shape analysis 

The spectra of many infrared-absorbing gases consist of 
individual lines whose shapes are pressure and temperature 
dependent. We have investigated the information contained in 
a single line observed with infinite resolving power. 

For an absorber in the lower part of the atmosphere the 
line shape can often be described by 

k (v) = , (7) 

tt ( Av + a (A) ) 

where S is the intensity, a is the collision-broadened width 
of the line, and k(v) is the absorption coefficient at 
Av = | v - v Q | . 

It is assumed, initially, that the intensity and width 
are temperature independent, that 
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C(A) = 0(0) , 


( 8 ) 


and that the data consist of absorptance values A(v) obtained 
from the observed spectral transmittances T(v), where 

A (v) = 1 - T (v) . (9) 

Provided A(v)<<1, 

A ( v n ) = Z Aj (v n ) , (10) 

j 

where 

A j ( v n) = m j k j p j 5ji j ' dD 

The spectral absorptances thus depend on the desired 
mixing-ratios weighted by the local air density, path length, 
and absorption coefficient. An example of such a spectral line, 
shown in Figure 65, is due to an absorber with a constant 
mixing-ratio . 

Near line center kj (v Q ) -S/fraj and thus, from Eqs. (7), 

(3) , and (11) -the spectral absorptance is independent of air 
density. Conversely, for Av>>a, the spectral absorptance is 
proportional to the square of the air density. It is these dif- 
ferences which might allow the inference of the height distribu- 
tion of an absorber from a single spectrum. 
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Transmittance 


A J' (cm' 1 ) 


Figure 65: A weak Lorentzian spectral line calculated for 

line intensity of tr, half-width of 0.1 cm - -'-, 
spectral value spacing of 0.004 cm - - 1 -, a constant 
mixing-ratio of 1.0, and a ground-level observer 


Local mixing-ratios were retrieved by linear, least 
-squares analysis of 125 equally spaced A(v) values, over the 
range 0<Av<5a(0), to which normally-distributed noise was added 
to give a SNR of 100 at Av=°°. 

Some results for N R =2, 4, and 5 are shown in Figure 66. 
The uncertainty increases more slowly with altitude than in 
Figure 18 but, at ground-level, it is an order of magnitude 
larger. Thus, even in this idealized situation, essentially 
no height discrimination is possible above the instrument. 

In practice, information is lost because of finite 
spectral resolution and the inverse exponential, rather than 
linear, change in absorptance with increasing absorber amount 
in the path. Also, in the mid-infrared, where many features of 
trace atmospheric constituents occur, the Lorentz width 
approaches the Doppler width at heights between 30 and 40 km. 
Above this height the dependence of the absorption coefficient 
at Av. = 0 on the air density, decreases to zero and there is 
little advantage over the measurement of the equivalent width 
of the feature. 

It can similarly be expected that the ability to discrimi- 
nate the height of the absorber, based on the temperature 
dependence of the line width, is minimal because of the slow 
variation of temperature in the lower atmosphere. Line inten- 
sities are also temperature dependent and Toth (5) has suggested 
that the atmospheric temperature profile can be obtained by 
observing CO 2 absorption lines arising from high energy states 


220 




U (Percent Uncertainty) 

Figure 66: The uncertainty in the mixing-ratio values 

retrieved from data shown in Figure 65 for value 
spacings of 25, 12, and 10 km. 
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in solar spectra. However, only a small fraction of the 
total number of observable lines of most molecules have a 
strong temperature dependence. Analyses based on these effects 
will tend to emphasize the lower, warmer regions of the 
atmosphere. It is unlikely that these temperature variations 
can be used to obtain good profiles. 


4.3 ANALYSIS OF ATMOS DATA 
4.3.1 General comments 

In the previous section we have shown that there is little 
added height discrimination for a ground observer by considera- 
tion of the line shape. This means that a mixing-ratio profile 
retrieved from a single line of a single spectrum obtained at 
ground-level can only contain a small number of discrete values, 
if no a priori assumptions about the shape are made. If 
several lines in the same spectrum are analyzed simultaneously 
the uncertainty in the retrieved mixing-ratio values should be 
reduced, but the vertical resolution is not expected to be 
improved. If, in addition, the spectrum analyzed has finite 
spectral resolution and if the spectral line parameters are 
uncertain the actual information content of the spectrum is 
less. Thus, from a single spectral line observed at ground- 
level it is expected that not more than one well-determined 
mixing-ratio value can be obtained. Therefore, unless the 
SNR is very high, for ground-level observations it is only 
possible to retrieve the total column amount of the absorber 
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(or one mean mixing-ratio value) ( 6 ) . The retrieval of 
one mixing-ratio value is equivalent to assuming that the 
profile to be retrieved is constant with altitude. Since the 
mixing-ratio profiles of N 2 O and CH 4 are essentially constant 
with altitude below 15 km (7) , this assumption is appropriate 
for the retrievals discussed in this section. 

In October 1982 in Lexington, Massachusetts, several 
solar spectra covering the region 600 to 4700 cm - -*- were taken 
in a ground-level test of the ATMOS interferometer shown .in 
Figure 2. The path difference of the interferometer was 
103.92 cm and the spectral resolution was 0.0096 cm"*-. The 
average signal values as a function of wavenumber, shown in 
Figure 67, were obtained from 16 spectra in the region 2400 to 
2600 cm - *-. The wavenumber spacing of these data was 0.0075 cm - *- 
and the observer altitude was 0.5 km above sea-level. The 
experimental conditions for these spectra are summarized in 
Table 25. 

Since, as noted previously, not more than one mixing-ratio 
value can be expected from analysis of a ground-level spectrum 
we have analyzed the spectral data of Figure 67 by assuming the 
absorber mixing-ratio values to be retrieved were constant with 
altitude. By this method we have been able to obtain mean 
mixing-ratio values for N 2 O and CH 4 using the spectral curve- 
fitting technique described below. 
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Figure 67: ATMOS spectra taken in October 1982 in Lexington, 

Massachusetts. 
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Figure 67 



Table 25: Experimental conditions for the recent test of 

the ATMOS instrument. 


Location 

Time 

Altitude 

Apparent zenith 
angle 

Spectral region 
covered 

Path difference 
of instrument 

Spectral resolution 

Spacing of spectral 
values 


Lexington, Massachusetts 

October 1982 

0.5 km above sea-level 

60° 

600-4700 cm" 1 

103.92 cm 

0.0096 cm -1 
0.0075 cm' 1 
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4.3.2 The spectral curve- fitting technique 

The observed signal I(vj_) obtained at frequency Vj_ in 
an absorption spectrum is usually related to the absorber 
amount u encountered along the homogeneous path and the 
absorption- coefficient k(Vj_) of the absorber by Beer's law 


I (vi) = I 0 (\>i)exp[-k(vi)u] , 


( 12 ) 


where I 0 (vi) is the signal obtained in the absence of the 
absorber and u is given by m x a where m is the mixing-ratio 
in the homogeneous path and a the optical airmass. 

In the lower portion of the atmosphere the line shape can 
be assumed to be Lorentzian and given by 


k i(vi) = 


a, 


< v i " ^o *) 2 + a ? 


(13) 


for the £th spectral line. The dependence of the half-width 
a on the temperature and pressure is often written as 




a 


o2. 




(14) 


and the line intensity as 


■S 


l ~ 



6 

exp [ -E^ (T 0 - T 


)/0 


6946T T 0 ] 


(15) 
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where 8=1 for N 2 O and CH 4 ( 8 ). The reference temperature 
and pressure T 0 and p Q for the AFGL line listing, used in this 
analysis, were 296 K and 1 atm respectively (9). 

In order to analyze the data in Figure 67 the path of 

the observed radiation was partitioned into three regions in 

which the atmosphere was approximated by a homogeneous layer. 

This was accomplished by dividing the ray path into three 

portions such that the optical airmass encountered in each of 

the portions was identical to the others . The regions thus 

selected were 0.5-3. 5 km, 3. 5-8. 5 km, and 8.5-40 km. The 

absorption coefficient was obtained assuming that the tempera 

ture in the ith region was given by the mass weighted average 
N i 


T-: = 


Z 

j=N oi 


a j T j 




(16) 


where N 0 j_ is the lower boundary of the ith region and the 


upper boundary, and the pressure was given by 


N. 

I 1 


Pi „ 3=N 


a j P j 


01 


(17) 


The optical airmass a^ in the ith region is 


N i 


a i = 2 a. 


j=N 


01 


where aj is the airmass in the jth layer. 


( 18 ) 
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It is convenient to define a coefficient K(v^) by the 
relation 


3 

K (v . ) = E k (v.)ar . (19) 

j=l 3 1 3 


Thus the signal for an absorber with constant mixing-ratio 
m is given by 


I(v i ) 


3 

I 0 ( v i) ex P[“ I ^( £ k. 

j=l 11 


(v i )ar)] . 


( 20 ) 


In general the absorption coefficient k^ is the sum of the 
contributions of all the lines in the spectrum. In practice 
an individual line will only affect- the spectrum over some 
wavenumber range about the line center v q . Beyond this range 
the contribution of the line will be below the level of the 
experimental noise. In the calculations for this analysis we 
have assumed a range of 10 cm ^ for all the lines in the 
spectrum. 

The spectrum analyzed, shown in Figure 67, contained two 
strong ^0 bands and some weak CH^ lines. We have attempted 
to retrieve both the ^0 and the CH^ mixing-ratio value assuming 
the mixing-ratio profiles to be constant with altitude. It was 
found that the spectral data were shifted by a constant amount 
over the entire spectral region. This shift was corrected 
prior to the retrieval of the mixing-ratio values. 
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Due to the strong wings of nearby CO 2 and Nj bands the 
background in the observed spectrum is not constant. If small 
wavenumber regions (say 10 cm are considered the background 
I 0 can be approximated by 

I o (v i ) = a o + a l (v i " v a) (21) 

where refers to the first wavenumber value in the 10 cm - '*' 
region. Thus the signal at v can be modeled by 

I (v i) = ( a o + ai(Vi- v A ))expr-m N 2 oK N2 o + m cH4 K CH4 )] 

( 22 ) 


if infinite spectral resolution is assumed. Lin et al (10) 
have shown that small frequency shifts in the spectrum can be 
modeled by including in the exponential the term 


A(v i ) = 


K (v i+l) - K(v i-1 ) 
v i+l~ v i-l 


Thus Eq. (22) becomes 


(23) 


KV = (a O + a l^i -y A ))eXp[ - (m N 2 0 K N 2 0 +m CH 4 K CH 4 


+ a o (A N 2 0 m N 2 0 + A CH 4 m CH 4 } } 1 


(24) 


230 



where a Q is the magnitude of the frequency shift. The 
retrieved magnitude of a 0 was always less than the wavenumber 
spacing of the data. 

At each wavenumber value the coefficients K. T _, , 

An 2 0 7 an< ^ ^ch were computed using the AFGL line listings (9) 
and assuming that the atmospheric temperature profile was 
given by the U.S. Standard Atmosphere (11) and the sea-level 
pressure was 1 atm. The parameters a Q , a^, a Q , m N Q , and m CH 
were obtained by fitting Eq. (24) to the spectral data shown 
in Figure 67 using non-linear least squares analysis techniques 


( 1 ) . 

The effects of the instrument, which are usually modeled 
by the convolution of Eq. (24) with an instrument function, 
were not included since the spectral resolution of the data 
shown in Figure 67 (the spectral resolution is given in Table 
25) was much less than the half-widths of N 2 0 or CH 4 at 
atmospheric pressure (9) . 


4.3.3 Some results from the analysis of ATMOS data 

For convenience the ATMOS data were divided into eight 
10 cm 1 regions as shown in Figure 68 . In each region the 
parameters in Eq. (24) which give the smallest RSS were obtained 
by the non-linear least squares analysis technique ( 1 ) . 

The residuals obtained for the retrievals performed on 
the eight spectral regions are shown in Figures 69-76. The 
largest features appearing in these residual plots are due to 
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Figure 68 
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Figure 69: The ATMOS spectra and the residuals obtained by 

spectral curve-fitting in the region 2430-2440 era' 
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Figure 70: The ATMOS spectra and the residuals obtained by 

spectral curve-fitting inthe region 2440-2450 cm' 
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Figure 71: The ATMOS spectra and the residuals obtained by 

spectral curve-fitting in the region 2450-2460 cm' 
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Figure 72: The ATMOS spectra and the residuals obtained by 

spectral curve-fitting in the region 2470-2480 cm' 
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Figure 73: The ATMOS spectra and the residuals obtained by 

spectral curve-fitting in the region 2480-2490 cm' 
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Figure 74 : The ATMOS spectra and the residuals obtained by 

spectral curve-fitting in the region 2490-2500 cm' 
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Figure 75: The ATMOS spectra and the residuals obtained by 

spectral curve-fitting in the region 2520-2530 cm 
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Figure 76: The ATMOS spectra and the residuals obtained by 

spectral curve-fitting in the region 2530-2540 cm' 



methane and water vapor lines not included in the AFGL listing 
and methane line intensities and positions with erroneous 
values. Another source of the features in these plots may be 
due to non-linearity of the instrument which causes the 
observed signal to increase non-linearly for increasing inci- 
dent radiation on the instrument. This effect was found to be 
more noticeable in the stronger lines in the spectrum (12) . 

Several retrievals, not shown, were performed on the 
portion of the spectrum above 2440 cm - "'" but the poor fit 
obtained in this region was attributed to the instrument non- 
linearity. 

The mixing-ratio values obtained for N 2 0 and CH 4 in the 
eight regions are listed in Table 26. Due to the poor values 
given by the AFGL listing for the methane line parameters, the 
mixing-ratio values obtained for CH 4 vary widely over the 
regions studied and the uncertainties of these values were 
large. However, consistent values were found for the N 2 0 
mixing-ratio values. The one exception was for region 7 
(2520-2530 cm' 1 )-. 

Listed in Table 27 are some values for the mixing-ratio 
value of N 2 0 obtained by other investigators. The retrieved 
mixing-ratio values listed in Table 25 were found to be in 
agreement with those found by other investigators, thus demon- 
strating that the non-linear, least squares technique is 
capable of extracting information about the absorber from the 
spectrum. 
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Table 26: Mixing-ratio values retrieved by spectral curve- 

fitting of the ATMOS spectra. 


Region 

# 

Av (cm - -'-) 

m N 2 0 x 10 “ 6 

m CH 4 x 10 " 6 

RSS 

1 

2430-2440 

0.257 ± 0.052 


0.091 

2 

2440-2450 

0.252 ± 0.026 


0.172 

3 

2450-2460 

0.271 ± 0.026 


0.109 

4 

2470-2480 

‘ 

0.266 ± 0.026 

1.30 ± 0.65 

0.272 

5 

2480-2490 

0.254 ±0.026 

0.45 ± 0.91 

0.086 

6 

2490-2500 

0.286 ± 0.13 

1.40 ± 0.52 

0.326 

7 

2520-2530 

0.187 ± 0.039 

1.22 ± 0.65 

0.411 

8 

2530-2540 

0.247 ± 0.026 

0.73 ± 0.39 

0.285 
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Table 27: Comparison of retrieved N 2 O mixing-ratio value to 

some values reported by other investigators. 


Retrieved 1 
of 

mixing-ratio 

N 2 O 

Experimental method 

262 

PPb 

average of retrieved values in 
this work (excluding region 7) 

264 

ppb 

value at 14.6 km from analysis 
of absorption spectra from 
balloon-borne instrument (13) 

30.8 

ppb 

air samples - aircraft (14) 

286 

ppb 

air samples - ground-level (15) 

255 

± 6 ppb 

air samples - aircraft (16) 

260 

ppb 

infrared absorption spectra 
from balloon-borne instrument 
(17) 

270 

PPb 

ground-level infrared absorp- 
tion spectral curve-fitting 
(18) 

310 

ppb 

ground-level infrared absorp- 
tion spectral curve-fitting 
(19) 

300 

± 10 ppb 

ground-level infrared absorp- 
tion spectral curve-fitting 
(20) 






Table 28: The effects of systematic errors on the mixing- 

ratio values retrieved from the ATMOS spectra for 
wavenumber region 2480-2490 cm - -'-. 


Physical 

parameter 

Systematic 

error 

Error in N 2 O 
mixing-ratio 

Error in CH 4 
mixing-ratio 

RSS 

surface 

pressure 

5% 

-2.0% 

-6.5% 

0.083 

temperature 

profile 

5 K at all 
altitudes 

2.0% 

2.2% 

0.088 

line 

intensities 

1% for all 
lines 

-1.2% 

-4.4% 

0.086 

line 

half-widths 

5% for all 
lines 

2.8% 

-2.2% 

0.083 




The effects of systematic errors on the retrieved 
mixing-ratio values were examined by analyzing the data in 
Figure 68 assuming that the atmospheric temperatures were 
higher at all altitudes by 5 K, that the sea-level pressure 
was 5% larger, that the line intensities given from the AFGL 
listings are in error by 1%, and the line half-widths in 
error by 5%. Each error was considered independently. The 
mixing-ratio values obtained for ^0 and CH^ together with 
the RSS retrieved for these systematic errors are listed in 
Table 28 for the spectral region 24 80-24 90 cm - '*'. Comparison 
of the uncertainties listed in Table 26 with those due to 
the systematic errors shows that the effects of these errors 
are less than the uncertainty due to the experimental noise. 

Thus it is expected that for spectra obtained at ground-level 
systematic errors should not be significant provided that the 
SNR<2 00 . 

These results show that in order to improve the retrieved 
mixing-ratio values from spectral data better spectral line 
parameters are needed in order to calculate the absorption 
coefficients. The residuals could also be reduced by including 
instrumental effects such as finite resolution and non- 
linearity in the retrieval calculation. The results of Table 
28 show that small uncertainties in the temperature and pressure 
values have little effect on the retrieved mixing-ratio values 
in this particular analysis. 
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4 . 4 SUMMARY 


In this chapter we have investigated the information con- 
tent of the shapes of spectral features. It was found that, 
for a ground-level observer, the height discrimination possible 
was not much better than that obtainable from the' analysis of 
equivalent width data and the uncertainty at ground-level was 

actually larger than that obtained from equivalent width 

/ 

analysis (for a layer spacing of 25 km and an SNR of 100 the 
uncertainty in the retrieved mixing-ratio at ground-level was 
12%) . This uncertainty is due in part to incorrect AFGL 
parameters (9) . Thus very limited vertical resolution is 
expected and only one mean mixing-ratio is expected from the 
curve-fitting of ground-level spectral data with an SNR of 100. 

Spectral data obtained from a recent ground-level test 
of the ATMOS ins trument were analyzed by non-linear least 
squares techniques in order to retrieve the mixing-ratio values 
of N 2 O and CH^. The N 2 O mixing-ratio value of 262 ppb obtained 
by this analysis was found to be essentially in agreement with 
mixing-ratio values reported by other investigators. 

For ground-level observations, systematic errors in the 
parameters describing the structure of the atmosphere or the 
spectral line parameters were found to have only a small effect 
on the retrieved mixing-ratio values (not more than a 3% error 
in the retrieved mixing-ratio value for a 5% error in sea-level 
pressure, 5 K error in the temperature profile, 1% error in the 
line intensities, or 5% error in the line half-widths) . The 
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uncertainty in the retrieved mixing-ratio value due to experi- 
mental noise was much larger (~10%) than that caused by the 
systematic errors, in agreement with the results of the previous 
chapters . 

The accuracy of the retrieved mixing-ratio values could 
be improved if instrumental effects such as non-linearity 
and the finite resolution of the instrument were included in 
the mathematical model. Improvement in these values (especially 
that obtained for CH4) could also be obtained if the spectral 
line parameters are more accurately determined and if the 
apparent zenith angle of observation is more accurately known. 

The analysis of ATMOS spectra obtained from ground-level 
is not a complete test of its possible performance in space. 

The data analyzed in this work were obtained from an early 
test of the ATMOS instrument which was conducted to determine 
any possible modifications needed in the experimental design. 
However, our analyses show that the ATMOS data are analyzable 
since the results in general confirm prediction. The spectral 
data collected by ATMOS in space are expected to be of much 
higher quality. Thus due to a more favorable observer geometry 
and higher quality data it is expected that accurate, high- 
resolution mixing-ratio profiles are retrievable from the 
analysis of ATMOS data. 
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Chapter V 
CONCLUSIONS 

5.1 INTRODUCTION 

In this work we investigated the retrieval of constituent 
mixing-ratio profiles from solar absorption spectra by con- 
sidering primarily the information contained in a single weak 
spectral line. This information consists of the total absorber 
amount, which can be obtained by measurement of the equivalent 
width of the line, and the shape of the spectral line. We 
investigated the re-trieval of mixing-ratio profiles from 
equivalent width analysis of spectral lines and curve-fitting. 

The methods used in this work to retrieve mixing-ratio 
profiles from spectral data required no a priori assumptions 
about the shape of the profile to be retrieved. These methods 
also retrieved all the mixing-ratio values simultaneously, 
and thus no errors were introduced into the solution by the 
analytical method. Our investigations have determined, for 
a given set of experimental conditions,' the limit to the 
quality of the retrieved mixing-ratio profiles. 

5.2 EQUIVALENT WIDTH DATA 

From the analyses of sets of 125 synthetic total absor- 
ber amounts, corresponding to 125 different paths through 
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the atmosphere and observer altitudes of 500, 30, and 0 km, 
calculated using the slant path program of Snider and 
Goldman (1) and a constant mixing-ratio profile, we 
retrieved up to 100 discrete constituent mixing-ratio values 
between 0 and 50 km altitude by the linear least squares 
technique. The program used for the retrieval (P9R (2)) 
obtained both the discrete mixing-ratio values and their 
standard deviations. The retrieved standard deviations were 
used to. determine the uncertainty in the mixing-ratio values 
due to simulated random noise as a function of altitude. 

The uncertainty in the retrieved mixing-ratio values was 
found to increase exponentially with altitude (see Figures 
16, 18, and 19), to be inversely proportional to the SNR 
(see Tables 4, 8, and 10) and Nq 35 (see Tables 5, 9, and 11), 
to be proportional to the number of retrieved mixing-ratio 
values N r (see Figures 16, 18, and 19), to be independent of 
the observer altitude (Figure 19), and was 45 times larger 
for soundings in the atmosphere above the observer than 
below (Figure 18) . These general relations were true 
independent of the mixing-ratio profile -examined (section 
2.7) . 

The vertical resolution of the retrieved mixing-ratio 
profiles for an uncertainty constant with altitude degraded 
exponentially with altitude (Eqs. (2.83) and (2.84)), was 
inversely proportional to the SNR and N 0 2 , was proportional 
to the number of retrieved mixing-ratio values, was independent 
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of the observer altitude, and was 45 times worse above the 
observer than below. These vertical resolution relations can 
be used to determine the spacing of the mixing-ratio values 
which gives an uncertainty in the retrieved values constant 
with altitude (section 2.6.2). 

The uncertainty expressions obtained in this work sug- 
gest that many investigators may have underestimated the 
errors in their reported retrievals from solar spectra. By 
the methods outlined in this work it is now possible to 
estimate the uncertainty in a mixing-ratio profile retrieval 
based entirely on the experimental conditions and not on any 
a priori assumptions about the errors. 


5.3 SYSTEMATIC ERRORS 

We have also investigated the effects of systematic 
errors in the sea-level pressure, the temperature profile, the 
observer altitude, the apparent zenith angle of arrival of the 
radiation, and the latitude of observation on the retrieval 
of mixing-ratio profiles from absorber amount data. 

A synthetic data set with no added noise and a constant 
mixing-ratio profile was analyzed assuming errors in each of 
these quantities. The systematic error as a function of alti- 
tude was obtained by comparison of the retrieved values to 
the true profile. The errors obtained agree with those pre- 
dicted from variation of the absorber amount equation. 
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The mixing-ratio error due to temperature profile error 
was found to increase linearly with altitude and was zero 
at an altitude approximately equal to the scale height 
(Figure 59) . The mixing-ratio errors due to errors in the 
sea-level pressure, the observer altitude, the apparent 
zenith angle, and the latitude were found to be essentially 
independent of altitude (Figures 58, 60-62) suggesting that 
the effects of these systematic errors on the mixing-ratio 
profile are indistinguishable. 

Other errors due to neglect of refraction, motion of the 
observing instrument and the field of view of the instrument 
were also investigated. The error due to neglect of 
refraction was found to increase- rapidly with decreasing alti 
tude below 25 km (Figure 57) . The smearing due to motion of 
the observer causes the optical airmass traversed to change 
by up to 30%/sec (Figure 51) and the smearing due to the 
field of view of the instrument implies that the observed 
radiation will have traversed an optical airmass uncertain by 
up to 70% (for a 2 milliradian field of view) . These effects 
must be properly modeled in order to retrieve meaningful 
mixing-ratio profiles. 

An internal calibration method was presented by which 
these smearing effects and the other systematic errors can 
be eliminated. The mean optical airmass traversed by the 
radiation (as well as the tangent height of the radiation) 
can be determined by measuring the equivalent width of a weak 
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spectral line due to a gas with known constant mixing-ratio 
such as CC >2 in the spectrum containing the constituent spectral 
line to be measured. Effective values can be determined for 
•the apparent zenith angle and the observer altitude from 
this optical airmass which eliminate the systematic errors. 

In addition, since both the C0 2 line and the constituent line 
are subjected to the same weighting over the path due to the 
smearing effects, these effects are eliminated. This allows 
many of the effects due to systematic errors and smearing to 
be removed from the retrieved mixing-ratio profiles of trace 
gases. 

5.4 SPECTRAL CURVE-FITTING 

The additional information contained in the shape of the 
spectral line was examined by using a non-linear least squares 
technique (2) to fit the spectral line shape. An ideal line 
was considered with no degradation due to finite resolution 
or smearing effects. It was found that for a Lorentz-shaped 
line the uncertainties in the retrieved mixing-ratio profile 
increased more slowly with altitude than that obtained from 
analysis of equivalent widths (Figures 18 and 66) but was 
much larger at 0 km for an SNR of 100. Thus the vertical 
resolution will not be much better if the spectral shape is 
analyzed than that if it is obtained from the analysis of 
equivalent widths. 
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A curve-fitting technique was used to analyze spectra 
near 2500 era - '*’ obtained from a recent ground-level test of 
the ATMOS instrument to obtain mean mixing-ratio values for 
N 2 O and CH^. The retrieved CH^ mixing-ratio value was not 
accurate due to poor line parameter values for CH^ in the 
region studied given by the AFGL line listings (3) . The 
mixing-ratio value obtained for N 2 O was in agreement with 
those reported by other investigators (Table 27) but uncer- 
tainties in the apparent zenith angle and possible uncorrected 
non-linearities in the data preclude obtaining definitive 
absolute values for the N 2 O mixing-ratio. 

5.5 FINAL COMMENTS 

The uncertainty and the vertical resolution of mixing- 
ratio profiles retrieved from spectral analysis depend on 
the number of observations, the tangent height spacing of 
these observations, the SNR, and the method of analysis. 

For an SNR of about 100 the vertical resolution of mixing- 
ratio profiles obtained from equivalent width analysis of 
a single spectral line from ATMOS spectra will not be better 
than 5 km due to spectral smearing. Thus the spacing of the 
retrieval values should not be less than 2.5 km. The 
uncertainties in the retrieved values will increase exponen- 
tially with altitude for a constant value spacing, however, 
these uncertainties can be reduced by increasing the mixing- 
ratio value spacing or by analyzing a number of spectral 
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lines in each spectrum. 

Finally C0 2 lines in each spectrum can be used to inter- 
nally calibrate the ATMOS data to eliminate the effects of 
the systematic errors and the spectral smearing. The resul- 
tant uncertainty in the values retrieved from equivalent 
width analysis of a single spectral line is not expected to 
be less than 10% due to uncertainties in the spectral line 
parameters (line intensities and half-widths) necessary for 
the internal calibration. However, this uncertainty can be 
reduced by using more accurate spectral parameters. 

Other analytical techniques such as the simultaneous 
spectral curve-fitting of several spectra together with 
more accurate spectral line parameters may allow mixing-ratio 
profiles of trace gases with high vertical resolution and 
small uncertainties to be retrieved from analysis of ATMOS 
spectra. The ATMOS experiment will allow a large vertical 
range of coverage, the simultaneous retrieval of the mixing- 
ratio profiles of many trace gases, and the continuous 
monitoring of the temporal variations of these profiles. 

This information is essential for an understanding of the 
atmosphere . 
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Appendix A 

THE THEORY OF SLANT PATH CALCULATIONS 
A. 1 INTRODUCTION 

To analyze spectral data obtained from radiation 
passing through the atmosphere it is necessary to know the 
influence of the atmosphere on the path traversed by the 
radiation and the amount of absorber encountered by the 
radiation. The variation of the air density and the 
refractive index of air with altitude make it difficult to 
do these calculations accurately. If the atmosphere is 
divided into concentric layers , as shown in Figures 3 and 
77 , the total path through the atmosphere is the sum of the 
paths through each layer. 

The following derivation for the path traversed and 
the optical airmass encountered in each layer is based on 
the analysis of Treve (1,2). These relations have been used 
in the program described by Snider (3) and Snider and 
Goldman (4) . A version of this program, modified to allow 
observer altitudes up to 500 km was used to perform all the 
ray path calculations in this work. 
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Figure 77: The slant path of radiation through the refracting 

atmosphere . 
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A. 2 ASSUMPTIONS USED IN THE DERIVATION 


The method described by Treve makes several assumptions 
about the physical properties of the atmosphere and the 
shape of the earth. The effects of turbulence in the 
atmosphere have not been considered. The assumptions used 
were : 

1) The local surface of the earth is spherical, 
although the earth's radius is latitude depen- 
dent. 

2) The index of refraction of the atmosphere 
depends only on the altitude. 

From assumption 1) the radial vectors joining the 
earth's center and the points along the ray path (shown in 
Figure 77) will be normal to the local earth's surface. 

Based on assumption 2) the path will lie in the plane 
defined by 0, J. and the center of the earth (see Figure 77). 
The ray tracing problem reduces to a two-dimension'al calcu- 
lation in this plane. These assumptions also imply that 
Fermat's principle yields 

nrsin z = constant (1) 

where z is the local zenith angle of the ray, r the radius 
vector, and n the index of refraction. 

3) The index of refraction is given by the 
Lorenz-Lorentz formula (5). 
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4) The atmosphere contains no water vapor. 

5) The composition of the atmosphere is constant 
with position so that the number of dispersive 
molecules per unit volume is proportional to 
the density. Thus the Lorenz-Lorentz formula 
is 


_1_ n 2 - 1 
P m n 2 + 2 


constant 


( 2 ) 


where p m is the mass density of air. 

6) The atmosphere is in hydrostatic equilibrium. 

Consequently the pressure P is given in terms of 
the altitude A by the differential equation 

dp = -gp m dA (3) 


where the acceleration due to gravity g is alti- 
tude dependent. 

7) The atmosphere is a perfect gas obeying 


PM 

Pm r*t * 


(4) 


8) The pressure is given by 

P = kpi" . 
m 


( 5 ) 
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By use of the other assumptions r is found to depend 
on the atmospheric temperature gradient which varies with 
height. Although Eq. (5) is used by Treve (1,2) and Snider 
and Goldman (4) it is not necessary for the derivation of 
the relations used for the slant path calculation, as will 
be shown. It is used here to be consistent with the 
formalism developed by Treve and Snider and Goldman. 

9) The variation of the index of refractions with 
wavelength of dry air at T = 288 K and P - 1 atm 
is given by Edlen's formula (6) 

(n-rl) xlO 8 = 6432.8 + 2949810. (146-1/X 2 ) -1 

2-1 ^ ® ) 

+ 25540 (41-1/X*) 

where the wavelength X is in ym. 

10) The ith layer is characterized by the values A^ 
and A^ + ^ of the altitude at the bottom and top of 
the layer, the temperature T^ at the bottom of 
the layer and (dT/dA) ^ within the layer. 

A. 3 THE GEOPOTENTIAL ALTITUDE 

It is convenient for calculations describing the earth's 
atmosphere to use the geopotential height equivalent for the 
geometric altitude A. The potential energy per unit mass at 
altitude A is given by 

y = /o g(A)dA . (7) 
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The geopotential meter m' is defined as the work done in 
lifting a unit mass one geometric meter through a region in 
which the acceleration due to gravity is identically 
g Q = 9.80665 m/s. Thus the geopotential height is given by 


U 

H = — - (8) 

9o 


2 2 

where g Q ' has magnitude g Q and units m /s m'. Combining 
Eqs. (7) and (8) gives 


H 



g(A)dA . 


(9) 


The acceleration due to gravity at altitude A is 
obtained using Newton's law of gravitation 


g 


9o 


Re 

Re+A 


2 


( 10 ) 


The geopotential equivalent of geometric altitude A is 
found by evaluation of the integral 


H 



\ Re+A | 


2 

dA . 


Integration of Eq. (11) yields 


(ID 


H 


go R e A 

g Q ** Rq+a 


( 12 ) 
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The term g Q /g 0 ' ( = 1 m^/m) is the conversion factor from 
geometric to geopotential meters. If y = g 0 /g Q ' Eq. (12) 
becomes 


H 


( R e A 

Y Ir^Ta 


( 13 ) 


In terms of the radial position r ( = R Q +A) the geopotential 
height is 



(14) 


Using Eq. (13) the geopotential height of any geometric 
altitude can be calculated. The zero of geopotential is at 
sea level at all latitudes. 


A. 4 THE RAY GEOMETRY 

The path of radiation leaving the source at astronomical 
zenith angle and arriving at the observer at apparent 
zenith angle Z Q is shown in Figure 77.* By similar triangles 
the angular refraction W is 

W = Z +n - Z Q (15) 

where Z is the zenith angle at intermediate point J. 

\ 

*Note that any point along the path can be considered the 
source or the observer, so that Figure 77 and the formalism 
developed here are also true in each individual layer. 
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If the ray is traced towards the source from point J 
a distance ds , as in Figure 78, the radius r is increased 
by dr and the polar angle by dfi . Thus 


tan Z 


rdft 

dr 


and 


(16) 


cos z 


dr 

ds * 


In terms of the zenith angle Z 


dfl 


tan Z 


dr 

r 


(17) 


(18) 


and 


ds 


dr 

cos Z 


(19) 


Application of the law of sines on the triangle shown 
in Figure 79 gives 


r o 

sin (Z 0 -9 ) 


r 

sin (180°-Z o ) 


( 20 ) 


where 6 is the polar angle of the apparent location of point 
J. Equation (20) simplifies to 
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Figure 78 ; 


Tile incremental 


siant path. 


267 




rsin(Z o -0) = r Q sin Z Q . (21) 

A. 5 THE ATMOSPHERIC PARAMETERS 

The hydrostatic equation can be expressed in terms of 
the geopotential height using the differential form of 
Eq. (9) 


gdA = g 0 'dH • 


( 22 ) 


Combining Eqs . (22) and (3) gives 


dP = -g 0 "p m dH . 


(23) 


The differential of Eq. (5) is 


dP = krp£' l dp m . 


(24) 


The ideal gas law together with Eq. (5) gives 


dT . Mkir - i)p r ~ 2 dp m 
R* 


(25) 


Substitution of Eqs. (22) and (23) into Eq . (25) yields 


g 0 M (r - l) 

dT = - T- -%-i^dH . 


(26) 
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Integration of Eq. (26) gives the temperature as 


T = T 0 + L q (H-H 0 ) 


( 27 ) 


where 


Lo 


AT 

AH 


g 0 "M ( r - 1 ) 
r* r 


(28) 


The exponent T of Eq. (5) is given in terms of the tempera- 
ture gradient by 


g 0 m/r* 

L 0 + g 0 " M / R * ' 


(29) 


Thus T varies with altitude and for the U.S. Standard 
Atmosphere 1976 (7) has values in the range -1.4<r<1.4. 

The expression given in Eq. (5) ,. which has the form of the 
adiabatic equation, does not correctly describe the earth's 
atmosphere since the atmospheric motion is generally not 
adiabatic (8) . As will be shown, assumption 8) is not 
necessary. Alternate derivations performed without using 
this assumption yielded identical relations for the pressure 
at altitude A. This derivation is described in the next 
section. The formalism using r is given below. 

Rearranging the terms in Eq. (24) gives 
dP 

dPfl, = vTp 1 " 1 (30) 

m 
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Combining with Eq. (23) , Eq. (30) becomes 


dp m 

~p7~ 


g c dH 


m 


kTp 


r-i 


m 


which with the ideal gas law reduces to 


dPm 

Pm 


1 9o' M dH 

r r* t * 


But T(H) has been found in Eq. (27). Let a 
8 = g 0 "M/R*T 0 . Equation (31) becomes 


dpm 

Pm 


If T ji 1 Eq 


0 dH 

r 1 + (H-H 0 ) 

To 

. (33) integrates to give 


ln = 6 ^olhU,+ ^- (H— H Q ) ] 

Pmo ■ " r L q T o 

or 


Pm 

Pmo 



(H-H c ) ]“ 


If T = 1 then L = 0. So Eq. (33) gives upon 


Pm = Pmo exp[- 3 (H-H 0 ) ] . 


(31) 

(32) 

: 1/(T - 1) and 

(33) 

(34) 

(35) 

integration 

(36) 
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With Eq. (28) and the definition of 8, Eqs 
can be written in the form 


9o M 
1 R*L„ 


P = { 


T q +L 0 (H-H 0 ) j 


P exp 
o c 


g 0 " M 


(H-H q ) 


Lq * 0 

Lq = 0 


(35) and (36) 


( 37 ) 


A. 6 ALTERNATE DERIVATION OF P(H) 

The hydrostatic equation, Eq. (23), when combined 
with ideal gas law gives 


dP- 

P 


g 0 " M 

R*T 


dH . 


(38) 


The earth's atmosphere up to an altitude of 86 km can be 
divided into segments in which the temperature gradient is 
constant. Thus the temperature is given by Eq. (27) . 
Substitution of T(H) into Eq . (38) gives 


dP _ M 9o dH 

P R* T 0 + L 0 (H-H 0 ) * 


If L 0 = 0 


P = P Q exp 


R*T 0 


(H-H q ) . 


(40) 


If L 0 ^ 0 then 


P = P, 


T 0 + L_ (H-H 0 ) 


»o M 

R*L n 


(41) 


272 



These equations are identical to Eq. (37). Therefore it is 
not necessary to assume the gas motion is adiabatic (assump- 
tion 8)), since the same relations are found by assuming 
that the temperature gradient is constant in discrete regions 
of the atmosphere. 


A. 7 THE POLAR ANGLE ft 

The polar angle Q of point J can be found in terms of 
the apparent polar angle 8 of point J from (see Eq. (16)) 


dJ2 = tan Z 


dr 


(42) 


where dr/r is given by the differential of Eq. (21) 


drsin(Z o -0) - rcos(Z o -0)d0 = 0 


(43) 


Therefore 


Y- = cot (Z 0 - 9)d0 


(44) 


Fermat's principle can be expressed as 


nrsin Z = n 0 r 0 sin Z Q . 


(45) 


273 



If this equation is combined with Eq. (21) 


n n 

sin 2 = — sin(Z o -0) 
n ° 


Thus 


tan Z = 


n D 

sin ( Z Q - 8) 


Vi- - 

v n 


sin^ (Z n - 0) 


Let 


u = 1 - 


n 

l n o / 


Equation (43) combined with Eq. (47) gives 


tan Z = cos (Z °~ 9) 


Vcos^ (z Q - e) - u 


tan (Z Q - 0 ) 


or 


tan Z = 


1 - 


u 


-1 


cos^ ( Z Q - 8 ) 


tan (Z 0 - 0 ) 


Combining Eqs . (^.2), (44), and (50) yields 


dQ = 


1 - 


u 


2 

cos ( Z Q - 0 ) 




de 


The apparent polar angle at the source 0 S is given from 
Eq. (21) by 


( 46 ) 


(47) 


(48) 


(49) 


(50) 


( 51 ) 
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( 52 ) 


. _i( R e + A o 

e s = Z 0 -S ln [ sin z 

where A Q is the observer height and A g the source height. 
Integration of Eq. (51) gives 

0 

n = f S 1 ! \ d0 . (53) 

J ° cos -(Z 0 -6)1 



A. 8 THE REFRACTIVE ANGLE W 

The angular refraction W is given in differential form 
from Eq. (15) by 

dW = dZ + dft (54) 

where dft is defined by Eq. (16). The differential of 
Eq. (45) gives 

nsin Z dr + rsin Z dn + nrcos Z dz = 0 . (55) 
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(57) 


dr 

r 


dn „ co t Z dz 
n 


Thus 


7 dn - dz 
dfl * ~ tan 2 n 


(58) 


and 


dn 


dW 


- -tan Z n 


(59) 


The 


differen 


tial 


o{ E q. («) g^es 


2 ndn 

du * - 2 * 

n„ 


Equations 


(48) and (60) y ie 


ield 


1 du 


{16 ), (59) r and (61) 

Combining Eqs. W' 


1 dft j: - du 
dW = 2 r dr i- u 


(60) 


( 61 ) 


(62) 


or 


1 r dfl du 
dW " 2 1 -u d© dr 


(63V 
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With Eqs. (20) and (51), Eq. (63) becomes 


dW 


1 r 0 sin Z Q 

2 ( 1— u) sin (Z Q - 8 ) 


u 


-h 


cos' 2 (Z n - 0) 


du 

dr 


d0 


where u and du/dr are functions of 8. 


A. 9 DETERMINATION OF u(_0) AND du ( 8) /dr 

The quantity u(0) and consequently du(8)/dr can be 
obtained using the Lorenz-Lorentz formula 

o _ 2 _ i 

1 _ n 2 - 1 _ _ 1 _ 

Pm n 2 + 2 p mo nQ 2 + 2 

and introducing the terms 


a = 1 


Pm 

p mo 


and 


c 

~0 


n. 


- 1 . 


Rearranging Eq. 

(1-u) - -ij 

n o 

2 

(1-u) - 2 

n o 

Solving Eq. (68) 


(65) with Eq. 



p mo no^ + 2 


for u 


(48) 


gives 


(64) 


(65) 


( 66 ) 


(67) 


( 68 ) 
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u 


(69) 


n 


1 2 
2 + 


= 1 + 


o 


Pm n o 


n p p mo n Q 2 + 2 


n, 


2 _ 


Pm 

p mo n Q 2 + 2 


- 1 


Combining the terms in Eq. (69) and multiplying by 
(n Q 2 + 2)/ (n Q 2 + 2) gives 


_ 0 0 n 2 + 2 

(n Q - 1) - (n Q + 2) + 


Pmo 


Pm (n f 2 


u = 


2 ~ 2 n v ° 

no p mo 


p m 

Pmo 


(n Q - 1) - (n Q + 2) 


Collecting terms in P m /p mo 


n. 


Pmo (n o 2 + 2 ) ( n o 2 -D “ (n 0 2 + 2) (n Q 2 - 1) 


Pm 

5 mo 


u = 


(n 0 2 -l) - (n 0 2 + 2) 


In terms of oo Eq. (71) is 


n, 


2 w(n 0 2 + 2) (n Q 2 - 1) 



u 

:= 







00 

(I - no ) - 3 


or 








00 

2 

(n 0 

- 1) (n Q 2 + 2) (n Q 

+ 1) 



n o 





u 


3 (1 

+ ^ co(n 0 2 - 1) ) 


Substituting 

£ o 

from Eq. (67) 

into 



2e 0 w 



2 



2 

• (1 

+ ~ 2 ~) (3 + 2e Q + £ 

o > 


u 

n o 





3 1 +|ue 0 (l + £ q /2) 


- 1 ) 


(70) 


. (71) 


(72) 


(73) 


(74) 
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or 


u = 


2e a Cl. + £o/2 ). Cl + j £ q Cl + e 0 /2 } ) / (1 + z 0 ] 2 
1 + j^£ q ^ + e o/2) 


(75) 


If 


= (1 + e, 


,) C1 + j 




-J -) ) / (1 


+ £ o ) 


(76) 


and 


c 2 


2 

"N 

J 



(77) 


then 


. u 


2e o(J j 


o^l + C 2 u) 


(78) 


The term e Q (the value of n Q - 1 at the observer) can 
be calculated using the Lorenz-Lorentz formula and the values 
of n and p m at sea-level. Thus e Q does not depend on A or 0 . 

The quantity du(0)/dr is found from the derivative of 
Eq. (78) 


du 

dr 


2C 1 C 2 £ o 
( 1 + (^u) 


du) 
7 dr 


But 


du 

dr 


d 

dr 



so 


( 79 ) 


2 79 



du _ 2 c 1 c 2 £ o d_ _Pm_ 

dr ’ (.1 + C 2 w) 2 dr Pm ° 


A. 10 T.HE OPTICAL AIPMASS 


The optical airmass in element ds of the ray path is 


da = p m ds 


in Kg/m . Combining Eqs . (18), (19), and (81) gives 


p m rda 

da = sin Z . 


Combining Eq. (82) with Eqs. (45) and (21) yields the result 


da = r Q sin Z Q p m 


(1-up d£ 
sin 2 ( Z 0 - 8 ) d0 


de . 


But pjjj can be expressed as a function of 0. The term p m (0) 
must be determined . before Eq. (83) can be evaluated. 


A. 11 DETERMINATION OF p m ( 6 ) AND ^-(Pm/Pmo 5 

The density equations (Eqs. (35) and (36)) can be 
simplified by defining 


H - H 0 
YR e - H o ' 
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From Eq. (14) h is given by 


r 



and using Eq. (21) 


h 


sin ( Z 0 - 8) 
sin Z Q 


Substituting in the variable h , Eqs . (35) and (3 

to 


Pm 

Pmo 


a -c 3 h) a r ji l 
e' c 4 h r = 1 


where 


c 3 


r - i 

r 


3(R e - H q ) 


and 


C 4 = 6 (Re - H q ) . 

The derivative S- is obtained from Eq. (32) 

dr Pmo 

d P m _ _ P m 1 9o M 1 dH 
dr p mo " p^o T R* T dr * 


(85) 

( 86 ) 

6) reduce 

(87) 

( 88 ) 

(89) 

(90) 


281 



Equations (84) and (85) give 


H - 

YR e - 



1 



Therefore 


dH 

dr 


^(yRe - H 0 ) 
r z 


or, using Eq. (21) 


(91) 


.(92) 


dH ( Y R e " H o) sin 2 ( Z Q - 8 ) 
dr = r o sin Z Q 


Thus Eq. (90) becomes 

■d Pm Pm 8 T o ( R ~ H o > sin2 (z o ~ 0 ) 

dr Pmo " Pmo F T r o sin Z Q 


A. 12 THE FINAL INTEGRAL FORMS 

The function u(0) can be determined by Eqs . (78), (66), 

(87), and (86) and du(9)/dr by Eqs. (80), (66), (87), (86), 

and (94) . Thus the differential expressions for dfi , dW and 
da can now be integrated. The polar angle between the source 
and observer is given by 





u(0) 


cos* (Z n - 0) J 


-h 


(95) 
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where 0 S is determined from Eq. (20) . The total angular 
refraction is given by 


w = 


F £ 


9 (Re - H n ) 


sin Z, 


1 - 


u 


-k T 


cos (Z Q - 6) 


o 

T 


(1-0)) sin (Z Q - 0) 

2 d0 . (96) 

(1 - u) (1 + C 2 oi) ^ 


And the optical airmass encountered in the ray path by 

j 

0 _ 

( 1 - 03 ) 


R sin Z Q 


a = 


mo e 
(R e - H 0 ) 


sin (Z D - 0 ) 


(1 - u) 



cos 



H 


d0 


(97) 


The expression given in Eq. (97) gives the optical air- 
mass in terms of the mass per unit area. The number of air 
molecules per unit area encountered along the ray is found 
by replacing the mass density p m by the number density p v 
which is given by 


P v = 


Pm N A 

M 


(98) 


where N ^ is Avagadro's number. The variable w in Eq. (97) 
will not change for this substitution since co is given 
in terms of a density ratio. 
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The number of air molecules per unit area in a vertical 


column is given by 
-kpop 

p v (r) dr . 


a v = 


R, 


(99) 


where R^p is the earth-centered radius of the top of the 
atmosphere. In terms of the geopotential height this 
expression can be written as 



( 100 ) 


If a v is defined as one airmass, then the number of air 
molecules encountered along the ray path can be expressed in 
terms of airmasses by taking the ratio a/a v (thus the term 
optical airmass) . 

Equations (95) -(97) were found by Treve to be well- 
behaved for altitudes less than 90 km. If the atmosphere 
is divided into a number of concentric layers, successive 
use of Eqs . (95) -(97) in each layer allows the total path 

of the radiation to be traced from observer to source. 

In each layer Z s becomes the angle of entry of the ray into 
the layer and Z Q the angle of exit. The optical airmasses 
calculated in this work were obtained using Eq. (97) . 

The theory presented in this appendix has been used to 
determine the temperature profile of the stratosphere (9) 
and the mixing-ratio profiles of NO, NO 2 t HNO^ , 0 3 , f^O, and 
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NC >2 (10-15) from the analysis of infrared solar absorption 
spectra. This formalism has been found to be useful in 
accurately approximating the paths traversed by solar 
radiation. 
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Appendix B 


Listed below are some of the symbols (and their 
definitions) used in this work. In general, subscripts . 

0, S, and T denote the value of the parameter at the observer, 
source, and tangent height respectively. The subscripts 

1, j, k, and l denote the value of the parameter in the ith, 
jth, kth, or Jlth layer. Exceptions to these rules are noted 


below. 

Symbol 

A 


Definition 

The altitude above sea-level in geometric 
meters. 


A v 

a 


a m 

a v 

F(Z q ) 


G 

g (a) 


The spectral absorptance at frequency v. 

The total optical airmass encountered 
along a ray path. 

The total amount of an absorber encountered 
along a ray path. 

The total airmass in a vertical column in 
the atmosphere. 

The transition energy corresponding to the 
£th spectral line. 

The fractional time rate of change of the 
optical airmass as a function of the apparent 
zenith angle. 

Newton's gravitational constant 
6.67xl0"^l nt - m^/Kg . 

The acceleration due to gravity as a func- 
tion of altitude. 
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Symbol 


Definition 


g Q The acceleration due to gravity at sea- 

level 

9.80665 m/s 2 . 

g Q ' 9.80665 m 2 /s 2 m^. 

H The altitude above sea-level in geopoten- 

tial meters. 


H 


s 


The pressure scale height in geometric 
meters . 


K(v) 


k 


v 


L 


i 





The effective spectral absorption coef- 
ficient for a layered inhomogeneous 
atmosphere as a function of frequency. 

The spectral absorption coefficient at 
frequency v. - 

The temperature gradient of the ith 
atmospheric layer. 

The path length of a ray through the ith 
layer of the atmosphere. 

The mean molecular weight of air 
28.966 gm/mole. 

The mass of the earth 
5.98 x 10 24 Kg. 


m (A) 


The mixing-ratio of an atmospheric con- 
stituent as a function of altitude. 


N The number of layers in the layered atmos- 

phere.. 

N a Avagadro's number 

6.02 x 10 22 molecules/mole. 


N R 


The number of retrieved mixing-ratio values. 


No 


The number of values in the retrieval data 
set . 


n 


i 


The index of refraction of air in the ith 
layer. 


n 


o 


The index of refraction at sea-level. 
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Symbol 
P (A) 

P(o) 

R e (4>) 

**o 

Res 


R* 


r 


s l 

ds 
T (A) 

Tv 

t 

*e 
U (A) 


u 

W 


w (am) 


Definition 

Pressure of the atmosphere as a function 
of altitude. 

The pressure at sea-level. Also used for 
sea-level pressure is P Q . 

The radius of the earth as a function of 
latitude . 

The vertical resolution of the retrieved 
mixing-ratio profile at ground-level. 

The vertical resolution of the retrieved 
mixing-ratio profile. Subscript A denotes 
the region above the observer and B the 
region below. 

The universal gas constant 
8.314 joule/mole K. 

The radial distance with respect to the 
earth's center of a point in the atmosphere. 

The spectral line intensity of the £th 
spectral line. 

Element of path through the atmosphere. 

Temperature of the atmosphere as a function 
of altitude. 

The spectral trasmittance at frequency v. 
Time . 

The period of circular satellite orbit. 

The uncertainty of retrieved mixing-ratio 
values as a function of altitude. 

1 - (n/n Q ) 2 . 

The refractive angle of a ray path bent 
through the atmosphere. 

The equivalent width of a spectral line 
corresponding to a ray path through absorber 
amount a m . 
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Definition 


Symbol 

Z 




A (Vi) 


£ 

0 

X 

U (A) 
v 

V A 


v 


oZ 

p v 


m 

a (A) 


y 

$ 


The astronomical zenith angle of a ray 
through the atmosphere. 

The apparent zenith angle of arrival at 
the observer of a ray through the atmos- 
phere. 

The half-width of the £th spectral line. 

The retrieved frequency shift obtained 
from spectral curve-fitting. 

The frequency shift coefficient of a Q at 
frequency Vi. 

Random experimental noise. 

Apparent polar angle of a point along the 
path of a ray through the atmosphere. 

The wavelength in jim of the observed 
radiation.. 

The gravitational potential energy per 
unit mass at altitude A. 

Frequency in cm -1 . 

The lowest frequency value of. the spectral 
data analyzed. 

Line center of the Zth spectral line. 

Volume density of the atmosphere. 

Mass density of the atmosphere. 

The standard deviation of the retrieved 
mixing-ratio values. 

The conversion factor from geometric to 
geopotential meters. 

Earth's latitude. 


n 


0) 


Astronomical polar angle of a point along 
the path of a ray through the atmosphere. 

1 - (p/p Q ) . 
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